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ABSTRACT 


Severe  difficulties  associated  with  the  design  of  satisfactory 
flush-type  HF  transmitting  antennas  for  supersonic  aircraft  have  prompted 
a  reconsideration  of  trai ling-wire  antennas  for  such  vehicles.  The 
aerodynamic  instabilities  which  caused  mechanical  failure  of  the  trailing- 
wire  type  of  antenna  and  led  to  its  abandonment  except  in  low-speed 
applications  are  carefully  examined  in  this  report. 

The  shedding  of  a  spiral  vortex  by  the  pear-shaped  and  spherical 
end-weights  customarily  used  was  primarily  responsible  for  rotary  motions 
which  occurred  at  short  cable  lengths.  This  behavior  can  be  eliminated 
through  the  use  of  a  streamlined,,  aerodynamics  1 ly  stable  end-weight. 

To  ensure  that  the  antenna  hangs  down  out  of  the  turbulent  boundary  layer 
and  wake  of  the  aircraft,  as  well  as  to  take  the  fullest  possible  advan¬ 
tage  of  the  favorable  electrical  properties  of  the  trailing  wire,  the 
end-weight  should  be  relatively  heavy  and  have  as  little  aerodynamic 
drag  as  possible, 

The  phenomenon  of  aerodynamic  wave  amplification,  which  leads  to 
violent  mechanical  cable  oscillations  in  the  vertical  plane,  is  subjected 
to  a  comprehensive  theoretical  analysis.  The  partial  differential 
equations  of  motion  of  the  curved  wire  are  shown  to  lead,  in  the  case  of 
harmonic  t ime -dependence ,  to  an  approximate  representation  in  terms  of 
three  waves.  A  general  method  is  given  for  solving  the  boundary  value 
problem  in  a  given  case.  Simple  analytic  expressions  are  obtained  for 
the  total  aerodynamic  amplification  of  a  traveling  wave  on  the  curved  wire. 
The  investigation  shows  that,  although  the  amplification  factor  cannot 
be  reduced  to  unity  at  supersonic  speeds  with  wires  of  known  materials 
and  weights  of  practical  sites,  it  can  nevertheless  be  held  below 
dangerous  levels. 

The  theory  is  applied  to  the  computation  of  certain  properties  of  an 
actual  experimental  trailing  wire  system  which  has  been  successfully  flown 
at  speeds  up  to  Mach  1.24.  An  amplification  factor  of  about  8  is  indicated 
for  the  worst  flight  condition. 
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Function  defined  by  Eq.  (B-?) 

Abbreviation  for  1  -  2i(y/v)(l  +  cr2)~^ 
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(1)  Gravitational  acceleration 

(2)  Amplitude  deviation  of  u 

Phase  deviation  of  u 

Non- aerodynamic  damping  constant 
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Aerodynamic  growth  or  damping  constant 
Mass  of  end-weight 

Dimensionless  transverse  displacement 

Dimensionless  transverse  displacement  (downstream  wave) 
Dimensionless  transverse  displacement  (upstream  wave) 
Dimensionless  transverse  displacement  at  particular  point 
Transverse  displacement 

Transverse  displacement  (downstream  wave) 

Transverse  displacement  (upstream  wave) 

Transverse  displacement  at  particular  point 
Maximum  transverse  displacement  in  standing  wave 
Transverse  input  amplitude  applied  at  point  £ 

Dimensionless  longitudinal  displacement 

(1)  Longitudinal  displacement 

(2)  Number  of  half-wave  loops  in  standing  wave 

Function  defined  by  Eq*  (B-41) 

Time-free  fractional  cable  tension  variation 

Total  cable  length 

Total  undistorted  cable  length 

Arc  length  along  cable 

Arc  length  along  undistorted  cable 

Cable  tension 
Equilibrium  value  of  T 

Time 

Wave  propagation  velocity 

Wave  propagation  velocity  in  vacuo 

(1  +  cr*  )K  w 

Homogenous  solution  for  u  (downstream  wave) 
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Dimensionless  vertical  coordinate  of  pofcnt  on  cable 

Variable  of  integration 

Mechanical  wavelength 

Mass  of  cable  per  unit  length 

Coefficient  of  viscous  friction  in  Af-direction 
Coefficient  of  viscous  friction  in  P-direction 
Dimensionless  frequency 

Dimensionless  natural  frequency  of  end-weight  and  suspension 
in  Af-direction 

Dimensionless  natural  frequency  of  end-weight  and  suspension 
in  P-direction 

Dimensionless  horizontal  coordinate  of  point  on  cable 
Air  density 

Total  dimensionless  cable  length 

Total  dimensionless  undistorted  cable  length 

Dimensionless  amplitude  of  second-order  length  variation 
Dimensionless  arc  length  along  cable 

Dimensionless  arc-length  coordinate  of  lower  end  of  cable 
Dimensionless  arc-length  coordinate  of  upper  end  of  cable 
Dimensionless  arc  length  along  undistorted  cable 

Dimensionless  time 
Slope  angle  of  cable 
Equilibrium  value  of  <f> 

Slope  angle  of  cable  at  lower  end 
Phase  angle 
Angular  frequency 


xii 


AERODYNAMIC  CHARACTERISTICS  OF  TRAILING. WIRE 
ANTENNAS  AT  SUPERSONIC  SPEEDS 


I  INTRODUCTION 


Most  aircraft  antennas  for  HF  transmission  act  merely  as  coupling 
devices  to  excite  currents  on  the  airframe;  the  actual  source  of  radiation 
is  the  airframe  currents.  This  is  true  of  tail  and  wing  cap  antennas,  and 
such  variants  of  these  as  wing,  tail,  and  nose  probe  antennas.1  It  is  also 
true  of  notch  and  shunt  antennas,  and  even  of  fixed  wires.2*3 

With  the  evolution  of  the  airplane  through  ever  more  sophisticated 
shapes  to  the  delta-wing  configuration,  however,  it  has  become  increasingly 
difficult  to  design  HF  antenna  systems  capable  of  achieving  adequate 
coupling.  While  on  the  one  hand  the  new  shapes  have  eliminated  many  of 
the  natural  charge  and  current  concentrations  necessary  to  obtain  effective 
coupling,  on  the  other  hand  the  weight  and  structural  penalties  which  can 
be  tolerated  have  become  ever  smaller.  Moreover,  the  drag  associated  with 
any  external  structure  has  risen  greatly.  The  problem  has  been  further 
complicated  by  the  fact  that,  in  the  few  remaining  possible  antenna 
locations,  high  ambient  temperatures  and  limited  space  make  extremely 
difficult,  if  not  impossible,  the  design,  housing,  and  maintenance  of  the 
often  quite  elaborate  antenna  couplers  which  are  required. 

These  difficulties  have  prompted  a  reexamination  of  one  of  the  older, 
more  primitive  types  of  antennas— the  trailing  wire.  This  antenna  is 
different  from  the  other  antennas  mentioned  in  that  the  wire  itself  is  the 
principal  radiating  element. 

Electrically,  the  trailing  wire  has  many  advantages  over  other  types 
of  HF  antennas.  One  of  the  most  important  of  these  is  that  it  can  be 
tuned  by  simply  varying  the  length.  The  magnitudes  and  degree  of  variation 
of  the  antenna  impedances  which  are  then  seen  at  the  feed  point  are  such 
that  very  simple  coupling  arrangements  will  suffice.  The  radiation  patterns 
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obtained  are,  in  general,  superior  to  those  of  the  other  antennas.  In 
particular,  if  the  end  of  the  wire  is  properly  weighted,  the  antenna  can 
be  made  to  hang  down  enough  to  give  considerable  vertical-dipole  radiation, 
thus  providing  more  nearly  onni -azimuthal  coverage.  In  addition,  because 
the  wire  itself  is  the  radiator,  both  patterns  and  impedance  tend  to  be 
independent  of  the  airframe  on  which  it  is  installed.  Although  the 
antenna,  when  extended,  contributes  appreciable  drag  to  the  airplane, 
it  can  be  retracted  when  not  in  use. 

In  spite  of  their  electrical  advantages,  trailing  wirec,  except  for 
special  low  speed  applications,  were  abandoned  some  years  ago;  at 
increased  speeds  aerodynamic  instabilities  set  in,  resulting  in  mechanical 
failure  of  the  wires  or,  in  some  cases,  damage  to  the  aircraft.  Since 
these  instabilities  were  responsible  for  the  abandonment  of  trailing  wires, 
it  is  appropriate  to  examine  possible  causes  for  them  to  see  whether  they 
might  not  be  eliminated.  Investigation  reveals  two  sources  of  instability 
which  appear  capable  of  accounting  for  the  observed  behavior. 

In  virtually  all  trailing-wire  antenna  installations  until  now, 
spherical  or  pear-shaped  end-weights  have  been  used.  It  is  known  from 
aerodynamic  investigations  that  a  sphere  in  a  high-speed  air  stream  sheds 
a  spiral  vortex.4  The  shedding  of  such  a  vortex  by  the  spherical  end- 
weight  on  a  trailing  wire  tends  to  induce  a  rotary  motion  of  the  sphere. 
Depending  on  the  length  and  tension  of  the  wire  a  standing  wave  can  be 
established  on  it,  resulting  in  large  amplitude  rotary  motions  of  the 
weight.  Such  motions  have  been  observed  in  actual  trailing  wires.  They 
occur  as  the  wire  is  reeled  in,  sometimes  causing  the  weight  to  pound 
against  the  bottom  of  the  aircraft  when  the  wire  becomes  short. 

Equally  serious  as  that  just  described  is  another  form  of  instability 
which  occurs  when  the  speed  of  the  aircraft  exceeds  the  speed  of  propaga¬ 
tion  of  mechanical  waves  on  the  wire.  When  this  condition  exists  a  wave 
traveling  downstream  on  the  wire  will  be  amplified  by  aerodynamic  action. 
This  phenomenon  was  investigated  by  W.  H.  Phillips  and  others.5'6  Phillips 
showed  that  the  downstream  wave  will  increase  exponentially  with  distance, 
while  the  upstream  wave  is  always  damped.  The  exponential  index,  which 
gives  the  amplification,  depends  in  a  complicated  way  upon  the  wind 
velocity,  the  amount  by  which  the  wind  velocity  exceeds  the  velocity  of 
propagation  on  the  wire,  the  angle  of  attack,  and  the  frequency  of  the 
disturbance. 
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Aerodynamic  amplif ication  of  the  type  jvst  described  accounts  very 
well  for  the  observed  failure  of  trailing -wire  antennas.  When  the  condition 
for  high  amplification  exists  a  snail  mechanical  motion  introdwced  near 
the  point  of  wire  attachment,  either  by  motion  of  the  aircraft  or  by 
turbulence  of  the  boundary  layer  air  through  which  the  antenna  passes, 
will  travel  downstream  on  the  wire,  arriving  at  the  end-weight  as  a  large, 
amplitude  ware.  At  the  end-weight  the  ware  is  partially  reflected  giving 
a  high  amplitude  standing  ware.  Moreorer,  coupling  through  wire  curratwre 
into  longitudinal  motion  leads  to  tension  rariations  in  the  wire.  The 
riolent  motion  existing  near  the  point  of  weight  attachment,  combined 
with  extreme  tension  rariations,  will  quickly  cause  failure  of  the  oable. 

The  rortex  shedding  responsible  for  instability  of  the  first  type 
described  abore  can  easily  be  dliminated  by  replacing  the  preriously  used 
sphere  with  a  streamlined  end-weight  deaigned  to  possess  inherent  aero¬ 
dynamic  stability.  Furthermore,  the  use  of  a  heary,  re ry  low-drag  end- 
weight  results  in  the  electrically  highly  farorable  cable  configuration 
in  which  there  is  a  region  of  appreciable  length  just  abore  the  end-weight 
where  the  wire  departs  only  slightly  from  the  rertical.  The  use  of 
stabilising  fins  does,  however,  give  rise  to  still  another  form  of  system 
instability  which  occurs  because  of  coupling  between  a  weather-rane  type 
of  motion  of  the  streamlined  end-weight,  and  a  lateral  pendulum  oscillation.7 
On  the  other  hand,  the  spatial  region  in  which  this  type  of  instability 
is  possible  is  so  close  to  the  aircraft  that  the  weight  will  never  normally 
be  in  it  except  when  passing  into  or  out  of  the  aircraft.  Furthermore, 
the  rate  of  growth  of  such  oscillations  is  so  small  that  they  cannot  attain 
significant  amplitudes  in  the  time  the  weight  spends  in  the  critical  region. 


Eliminating  the  second  type  of  instability,  that  due  to  aerodynamic 
amplification,  is  much  more  difficult.  Since  this  form  of  amplification 
occurs  only  when  the  component  of  wind  velocity  tangential  to  the  wire 
exceeds  the  velocity  of  ware  propagation  along  it,  our  first  thought  is 
to  see  whether  the  propagation  velocity  cannot  be  increased  to  a  value 
greater  than  the  wind  velocity,  even  at  very  high  speeds. 


To  a  first  approximation  we  can  say  that  the  propagation  velocity  of 
wares  on  the  wire  is  the  same  as  that  which  would  obtain  in  a  vacuum.  It 
can  be  shown  that  this  velocity,  t/0 ,  is  given  by 
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stress 

density 
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Hard-drawn  steel  wire  has  the  highest  available  ratio  of  ultimate  strength 
to  density  and  for  such  wire  stressed  to  the  yield  point  the  propagation 
velocity  of  the  wave  on  the  wire  is  approximately  Mach  2.  Since  in  any 
actual  installation  we  must  allow  a  substantial  factor  of  safety,  the 
highest  propagation  velocity  that  can  be  obtained  in  practical  antennas 
will  be  of  the  order  of  Mach  1.  If  the  airspeed  exceeds  this  figure,  the 
amplification  can  be  reduced  to  unity  only  by  increasing  the  weight  of 
the  terminating  body  to  such  a  value  that  the  entire  cable  is  nearly 
vertical.  For  aircraft  in  the  Mach  2  to  Mach  3  speed  range,  it  can  be 
shown  that  the  weights  theoretically  required  would  be  of  the  order  of 
10  or  20  tons,  and  the  cables  would  have  to  be  of  the  order  of  X  to 
1  inch  in  diameter.  Ve  therefore  conclude  that  it  is  impossible  practi¬ 
cally  to  avoid  having  some  amplification  when  the  trailing  wire  is  mounted 
on  such  aircraft.  Attention,  then,  must  be  directed  toward  keeping  this 
amplification  within  tolerable  limits. 

It  is  worthy  of  note  that  the  amplification  can  also  theoretically 
be  reduced  to  a  very  low  value  by  using  a  light  end-weight  of  very  high 
drag,  so  that  the  cable  streams  out  almost  horizontally.  In  practice, 
however,  this  design  places  the  entire  cable  in  the  turbulent  boundary 
layer  and  wake  of  the  air  vehicle,  where  the  airflow  is  certainly  not 
laminar  and  parallel  to  the  cable  at  all  points.  This  condition, , which 
may  very  well  result  in  extremely  high  amplification,  will  not  be 
considered  in  this  report. 

Phillips*  first  order  theory  applies  only  to  straight  cables  in 
uniform  airstreams,  and  takes  no  account  of  the  true  shape  of  the  antenna. 
The  actual  curved  shape  results  in  substantial  departures  from  first-order 
predictions,  particularly  at  lower  frequencies.  Furthermore,  Phillips* 
theory  assumes  constant  wire  tension.  Consequently  it  not  only  cannot 
predict  the  conditions  under  which  failure  occurs,  but  provides  no  insight 
into  the  circumstances  under  which  the  theory  ceases  to  be  approximately 
valid. 

Since  the  preceding  discussion  makes  it  clear  that  aerodynamic  wave 
amplification  cannot  possibly  be  avoided  at  speeds  significantly  above 
Mach  1,  it  is  important  to  have  a  theory  which  more  accurately  describes 
the  behavior  of  trailing  wires  at  velocities  where  amplification  can  occur. 
To  be  useful  the  theory  should  include  at  least  a  first  order  description 
of  the  tension  variations  which  can  occur  because  of  aerodynamic  wave 
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amplification,  aince  it  ia  obvious  that  wire  failure  cannot  take  place 
without  tension  fluctuations.  It  should  also  show  the  effect  of  the 
curved  shape  of  the  wire  upon  the  amplification  and  resultant  wire 
tension.  In  this  report  such  a  theory  is  developed  and  applied  to  an 
experimental  installation  designed  for  flight  testing  under  Air  Force 
Contract  AF  331616 ) -5549.  While  the  flight  tests,  because  of  limitations 
of  the  F-100  test  vehicle,  have  been  conducted  at  rather  low  supersonic 
speeds  (up  to  Mach  1.24),  the  installation  is  designed  to  simulate  the 
conditions  of  wave  amplification  likely  to  be  encountered  at  the  Mach  2 
to  Mach  3  speeds  representative  of  the  next  generation  of  aircraft. 
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II  ONE-DIMENSIONAL  APPROXIMATION  FOR  WAVE  AMPLIFICATION 


A.  PRELIMINARY  CONSIDERATIONS 

Before  investigating  the  dynamic  situation  in  which  wave  motion  exists 
on  the  trailing  wire  it  is  necessary  to  explore  the  static  condition, 
where  the  wire  is  in  equilibrium  under  the  action  of  the  aerodynamic  forces 
of  the  airstream  and  gravitational  forces  —  principally  those  due  to  the 
end-weight.  The  shape  of  a  trailing  wire  in  an  airstream  has  been  con¬ 
sidered  by  earlier  investigators.8,9  It  is  re-examined  here,  however, 
because  it  is  possible  to  derive  the  results  somewhat  more  succinctly  than 
has  been  done  heretofore,  and  in  a  form  more  suited  to  the  analysis  of  the 
dynamic  condition  which  follows. 


B.  DIMENSIONLESS  EXPRESSIONS  FOR  EQUILIBRIUM  WIRE  SHAPE 


4F 


T+dT 


The  aerodynamic  force  on  an  element  of  cable  suspended  in  a  supersonic 
airstream  is  always  nearly  normal  to  the  axis*  of  the  element,  and  can  be 
considered  to  arise  entirely  from  the 
component  of  air  velocity  in  the  normal 
direction,  V  .  10  For  cases  of  practical  in¬ 
terest,  the  weight  of  the  element  can  be 
neglected,  and  in  the  absence  of  mechanical 
waves  on  the  cable  the  element  is  in  equi¬ 
librium  under  the  action  of  the  aerodynamic 
force  and  the  tensions  at  its  ends  (see 
Fig.  1).  The  tensions  are  equal  in  the 
static  case  because  of  the  absence  of  any 
other  tangential  components  of  force — i.e., 
dT  «  0. 


FIG.  1 

FORCES  ON  AN  ELEMENT  OF  WIRE 


Let  the  following  quantities  be  defined: 
T  *  tension  in  cable 


*  TM  saall  coapoaeet  of  aerodyaaaic  fore*  taafeatial  to  the  cable  eleseat  baa  beta  aaflectad,  as  its  coa* 
side ratios  voald  bare  sada  the  aolutioa  of  aa  already  vary  coaplax  problaa  area  sore  iavolvad.  Tbe  pria* 
cipal  affect  of  tbia  force  is  to  cease  a  rariatioa  ia  tbe  static  teaaioa  a loaf  tbe  cable.  Sifaificaat 
departures  froai  the  value  predicted  by  the  siapler  theory  occur,  however,  oaly  ia  tbe  upper  part  of  loaf 
wires.  Ia  this  refioa  tbe  effect  of  iacreaeed  teaaioa  oa  cable  shape  ia  eii|ht.  A  secoad  effect  of  tbe 
iaervaaed  teaaioa— »a  decrease  ia  aerodyaaaic  asplificetioa— ia  also  quite  ssall. 


ds  ■  length  of  element 
4>  m  angle  element  makes  with  horizontal 
«  air  velocity  (horizontal). 

Then,  from  Fig.  1  one  can  write  for  the  mechanical  force  acting  on  the 
wire  element 

dF  -  -T  d4>  (1) 

The  aerodynamic  force  on  the  element  is 


pV\  PV\ 

dF  *  CD  — —  ( D  ds)  *  CqD  -  sin2  4>  ds  (2) 

2  2 

where 

Cp  *  drag  coefficient  of  a  cylinder 
D  *  diameter  of  cable 
p  •  air  density. 

For  equilibrium  the  mechanical  force  equals  the  aerodynamic  force,  so  that 


or 


-T  d4> 


sin2  <t>  ds 


ds 

T 


-esc2  4>  d<f> 


where  C  is  a  characteristic  length  given  by 


\pv\ 

I  2 


(3) 


(4) 


(5) 


Provided  that  the  drag  coefficient  CD  can  be  considered  constant  along  the 
cable, *  Eq.  (4)  may  easily  be  integrated  to  give 


*  S*«  Ste.  JV*A  for  i  diiciaaioa  of  tfco  iaportaica  of  tki*  liaitotioo. 
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cot  0  *  cot  •  +  a 


(6) 


as  the  equation  of  the  cable  shape,  where  a  is  defined  aa  the  dimensionless 
arc  length  s/C.  Although  not  apparent  when  written  in  this  form,  Eq.  (6) 
is  the  equation  of  a  catenary. 

The  slope  angle  of  the  lower  end  of  the  cable,  ia  determined  by 
the  way  in  which  the  cable  is  terminated  at  that  point.  If  the  termination 
is  provided  by  a  body  of  very  high  weight- to- drag  ratio  the  lower  end  of 
the  wire  will  be  essentially  vertical,  and  the  tension  will  be  equal  to 
the  weight  of  the  terminating  body.  If,  in  addition,  the  lower  end  is 
taken  as  origin  cf  dimensionless  Cartesian  coordinates  g  and  r),  defined 
as  x/C  and  y/C  respectively,  the  following  expressions  result: 


i 


da '  cos  0'  *  esc  0-1 


(7) 


V 


4> 

da'  sin  0'  ■  In  (esc  0  +  cot  0) 


(8) 


The  four  quantities  £,  rjt  a,  and  0  can  be  expressed  in  terms  of  each 
other  by  means  of  various  substi tutions  among  the  above  equations.  Com¬ 
plete  results  appear  in  Appendix  A,  and  curves  of  the  wire  shape  are  shown 
in  Fig.  2. 

C.  WAVE  AMPLIFICATION  IN  A  STRAIGHT  WIRE 

According  to  Phillips1  theory,  which  is  restricted  to  vibrations 
taking  place  in  the  plane  determined  by  the  air  velocity  and  the  wire 
direction,  the  amplification  depends  on  the  angle  between  wind  and  wire, 
and  on  the  ratio  of  wave  velocity  to  wind  velocity.  To  the  first  order, 
the  only  effect  of  the  wind  on  in-plane  waves  (for  sufficiently  high  fre¬ 
quencies)  is  to  cause  the  wave  amplitude  to  vary  exponentially  along  the 
wire.  For  a  wave  traveling  downwind  (negative-*  direction  along  the  wire 
in  Fig.  2),  the  transverse  displacement  for  a  sinusoidal  input  is  given  by 

n  .  a0  €Jk9  €**(•♦•/*)  (9) 
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where  V  is  the  were  Telocity  on  the  wire,  o0  is  a  constant,  and 

k  •  (1  — - — -Jain  20  .  (10) 

2 C  \  V0  cos  4>) 

It  is  seen  that  a  ware  traveling  in  the  negative-s  direction  will  be 
amplified  if  the  component  of  wind  velocity  along  the  direction  of  the 
wire  is  greater  than  the  wave  velocity  ( k  positive),  and  damped  if  it  is 
less  (k  negative).  For  a  wave  traveling  upwind,  only  the  sign  of  U  is 
changed  in  Eqa.  (9)  and  (10).  In  this  case,  k  is  always  positive  and, 
since  the  wave  travels  in  the  positive-*  direction,  damping  always  occurs. 

The  above  equation  indicates  that  there  can  be  no  amplification  when 
the  wire  is  exactly  parallel  to  the  wind  direction.  Likewise,  the  theory 
predicts  that  no  amplification  will  occur  for  vibrations  taking  place  in 
a  direction  normal  to  the  plane  containing  the  cable.  Actually  this 
result  stems  from  the  neglect  of  certain  terms  in  the  first  order  theory. 

A  more  accurate  investigation  of  wire  behavior  for  zero  angle  of  attack 
(see  Appendix  B)  shows  that  there  is  a  very  small  aerodynamic  amplification 
for  a  downwind  wave  and  very  small  damping  for  an  upwind  wave,  in  the  case 
where  one  or  the  other  exists  on  the  wire  alone. 

When  complete  reflection  of  the  lateral  wave  occurs  at  the  lower  end 
of  the  cable,  the  resulting  standing  wave  consists  of  a  combination  of 
the  two  traveling  waves,  but  the  net  damping  cannot  be  computed  by  means 
of  linear  superposition.  An  approximate  treatment  of  the  aerodynamic 
damping  is  still  possible,  however,  and  shows  that,  in  the  absence  of  other 
damping  mechanisms,  the  wire  behaves  as  a  very  high-Q  system  with  resonances 
at  multiples  of  one-half  wavelength. 

The  most  extreme  response  occurs  at  one-half- wave  resonance,  where 
it  is  found  that  the  standing  wave  amplitude  is  proportional  to  the  square 
root  of  the  input  amplitude  at  the  upper  end  of  the  cable.  The  theory 
developed  in  Appendix  B  takes  no  account  of  wire  curvature  or  of  other 
mechanisms  providing  coupling  to  the  in-plane  vibrations,  which  are  of 
much  lower  Q.  Moreover,  the  end-weight  does  not  constitute  a  perfectly 
rigid  termination,  and  therefore  will  furnish  some  damping  due  to  its 
lateral  drag. 

The  actual  Q  of  the  lateral  vibrations  will  therefore  be  considerably 
lower  than  that  predicted  by  Appendix  B.  For  this  reason,  and  more 
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importantly  because  first-order  coupling  to  tension  variations  occurs  only 
▼is  the  in-plane  notion,  it  ia  not  felt  that  the  lateral  motion  is  a 
significant  cause  of  cable  failure.  The  analysis  given  in  the  main  body 
of  this  report  will  therefore  be  confined  to  vibrations  in  the  plane  of 
the  wire. 

D.  TOTAL  AMPLIFICATION  IN  THE  CURVED  WIRE 

Equation  (9)  may  be  regarded  as  an  equation  for  a,  the  wave  amplitude 
along  the  wire,  which  may  in  the  general  case  be  complex: 

a  •  a9  e~k 9  (11) 

The  change  in  amplitude  in  an  element  of  the  wire  is 

da  m  - k  a#  €"**  ds  ■  -a  (cos  <p  “  y)  sin  <t>  dcx  (12) 

where 
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(13) 


An  approximation  to  the  total  amplification  in  a  curved  wire  can  be 
obtained  by  taking  Eq.  (12)  as  the  differential  equation  for  the  amplitude, 
and  inserting  the  known  relation  between  cr  and  <p  for  the  Equilibrium  shape 
of  the  wire.  This  is  equivalent  to  the  assumption  that  the  transverse 
wave  in  the  wire  travels  along  the  curve  as  if  the  wire  were  straight, 
remaining  always  transverse,  but  that  the  amplification  coefficient  varies 
in  the  correct  fashion  from  point  to  point. 

Since  dar  !.«  .tosc?  4>  d<fi,  we  have 

O'  i  ■  _i 

a*1  da  ■  (cos  <p  -  y)  esc  <t>  dcf>  .  (12) 


If  the  amplitude  at  the  origin  (lower  end  of  the  cable)  is  taken  to  be 
o# ,  and  U  (or  y)  is  assumed  to  be  independent  of  4>t  Eq.  (121)  leads  to 


(esc  4>  +  cot  4>)y  iin  4> 

(14a) 

£or  +  /l  +  a*Y  ( 1  +  cr2 ) 

(14b) 
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The  amplification  from  a  point  <r  on  the  wire  down  to  the  end  ia 


•  -  -  (1  +  <r2)^  [cr  +  /l  +  <r2] 

a(<r)  i 


(15) 


while  the  amplification  from  a  point  down  to  a  lower  point  •O’j  ia 

«(«*!>  A(  0,^) 


A(<rl,<r2)  ■ 


<*(<^2  )  dtO.'O'j  ) 


(16) 


E.  THE  REFLECTED  WAVE 

In  general  there  will  alao  be  an  upwind,  or  reflected,  wave.  The 
attenuation  experienced  by  such  a  wave  traveling  from  the  origin  to  a 
point  <r  can  be  found  from  Eq.  (15),  using  a  negative  value  of  y.  Loci 
constant  A(0,<7)  in  the  y,<r  plane  are  shown  in  Fig.  3.  For  example,  if 


FIG.  3 

AERODYNAMIC  AMPLIFICATION  CHARACTERISTICS  OF  TRAILING  WIRE 


of 
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the  wind  Telocity  is  twice  the  ware  Telocity  in  the  wire  (y  *  0.5),  the 
incident  waTe  will  be  amplified  by  a  factor  of  2.3  in  proceeding  fro* 

0  ■  10  to  the  end  of  the  wire,  while  the  reflected  ware  will  be  attenuated 
by  a  factor  of  45  in  returning  to  the  aa*e  point.  In  using  Fig.  3  it 
should  be  remenbered  that  y  •  0  corresponds  either  to  zero  ware  Telocity 
or  to  infinite  wind  Telocity. 

F.  THE  COMPLETE  SOLUTION 

The  conplete  solution  for  the  transTerse  displacement  on  the  cable 
for  the  case  of  a  single  sinusoidal  excitation  is,  in  the  one- dimensional 
approximation  of  this  chapter, 

al  ai 

n  9  -  ^  iu  (  t+Cc /  V)  +  -  €  io>  (  t-Ca/ V)  (17) 

Al(0,<r)  A210,<t) 

The  constants  a1  and  represent  the  (complex)  amplitudes  of  the  incident 
and  reflected  wares  respectively,  at  the  origin.  In  a/ given  case,  they 
can  only  be  evaluated  through  a  knowledge  of  the  appropriate  boundary 
conditions — on  the  one  hand  the  way  in  which  the  wire  is  terminated  at  its 
lower  end,  on  the  other  the  location  of  the  input  excitation  and  the  phase 

relationship  between  the  two  waves  at  that  same  point. 

! 

For  example,  if  no  transverse  motion  of  the  lowei  end  of  the  wire  is 
permitted,  then  a2  •  -Oj — i.e.,  complete  reflection  with  phase  reversal 
occurs.  If  the  input,  of  amplitude  nj,  is  applied  at  the  point  cr  ■  Z, 

that  is,  if  n(Z)  «  nj*  ,  Uj  can  be  evaluated  as  follows: 

ai  *  nl  t'4!1  *  Atl  e~iu>cl/v]--  1  .  (18) 

Vith  the  appropriate  substitutions,  Eq.  (18)  can  be  rewritten  as 

ttl  *  ^  n2  ^  +  2*  csch  sinh-1  2  +  “ — ^  (19) 

The  maximum  response  for  excitation  at  a  given  point  occurs  if  the 
frequency  is  such  that  that  point  is  a  node  in  the  standing  wave — i.e., 
the  distance  from  the  end  of  the  wire  is  an  even  number  of  quarter  wave¬ 
lengths.  If  this  distance  is  an  odd  number  of  quarter  wavelengths, 
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corresponding  to  excitation  nt  an  ant inode,  the  raaponaa  is  nininua. 

Used lass  to  say,  sines  the  reflected  ware  is  attennated  in  any  given  length 
of  wire  sore  than  the  incident  ware  is  amplified  in  the  ease  length,  the 
two  wares  can  nerer  cancel  each  other  at  aore  than  one  point,  in  this  case 
the  origin.  This  fact  tends  to  reduce  the  Magnitude  of  the  resonances. 

These  nay,  howerer,  be  quite  pronounced  for  snail  cr  (short  wires)  at  rela- 
tirely  high  frequencies  or  for  rery  snail  y  (rery  high  airspeed)  at  any 
frequency. 

The  foregoing  analysis  is  not  entirely  satisfactory  for  sereral 
reasons.  First,  the  effect  of  wire  shape  on  wares  whose  warelengths  are 
coaparable  with  the  radius  of  currature  has  not  been  taken  into  account 
in  any  way.  Second,  the  possibility  of  longitudinal  notion  of  the  wire, 
arising  fron  the  currature,  is  not  considered.  Third,  non-linear  effects 
of  rarious  kinds  are  entirely  ignored. 

G.  NON- LINEAR  EFFECTS 

The  non-linear  effects  Mentioned  abore  cone  into  play  when  ware  ampli¬ 
tudes  on  the  wire  are  not  entirely  negligible  in  conparison  with  warelength. 
One  such  effect  arises  fron  the  deriation  of  the  sine  function  fron  the 
tangent  function,  by  neans  of  which  it  is  approxinated  in  the  first-order 
theory. 

A  second  effect,  which  nay  hare  a  nore  profound  influence  on  the 
beharior  of  the  wire,  is  the  rariation  in  orer-all  length  caused  by  the 
presence  of  a  standing  ware.  Note  that  the  total  ware  on  the  wire  is  a 
conbination  of  a  standing  and  a  trareling  ware— alnost  conpletely  standing 
at  the  lower  end  for  a  rigid  ternination  there,  and  predoninant ly  trareling 
at  the  upper  end  in  noat  caaea  of  physical  interest. 

The  longitudinal  notion  at  the  lower  end  produced  in  this  way  will 
gire  rise  to  tension  fluctuations  if  it  is  necessary  to  accelerate  a 
relatively  large  nass.  It  should  be  noted  also  that  these. fluctuations 
will  take  place  at  double  the  ware  frequency.  If  their  anplitude  is  greater 
than  T,  the  cable  tension,  the  cable  will  go  slack  during  a  portion  of  the 
cycle,  then  go  taut  again  quite  suddenly.  In  an  extrene  case  the  cable  is 
slack  noat  of  the  tine  and  has  very  high  tension  for  an  extrenely  short 
interval.  This  condition  ia  obviously  one  possible  cause  of  cable  failure. 

It  ia  shown  in  Appendix  C  that  the  anplitude  of  the  length  rariation 
produced  by  a  sinusoidal  input  on  a  wire  of  length  £0,  when  the  incident 
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wave  haa  amplitude  ax  at  the  lover  end  of  the  wire  and  ia  completely  re* 
fleeted  there,  is  given  approximately  by  the  following  equation; 

Z'  »  20r|al|/x)*  tan"1  I0  *  (C-ll) 

In  the  above  equation  a}  and  k ,  the  wavelength  on  the  wire,  must  be  speci 
fied  in  terms  of  the  same  units,  and  £'  is  then  given  in  the  usual  dimen¬ 
sionless  units  of  this  theory.  Eq.  (C-ll)  holds  closely  only  for 
wavelengths  short  compared  with  the  radius  of  curvature  oH  the  wire. 

It  is  not  difficult  to  show  that,  if  the  input  consists  of  several 
sinusoids  having  frequencies  a>m ,  phase  angles  and  amplitudes  a1(|  at 

<J  »  0,  the  length  variation  may  be  expressed  by 


tan 


-1 


AI 


2'  cos  (2 co  t  +  2^_) 


or 


Q)* 


L  u 2 


cos  (20)nt  +  2 \pm) 


(20) 


An  extension  of  Eq.  (20)  leads  to  the  interesting  corollary  that  the 
frequency  spectrum  of  the  input  must  have  delta-function  type  singularities 
in  order  that  there  be  any  effect.  A  “white”  input  produces  no  variation 
in  length.  This  fact  is  a  result  of  the  dependence  of  Eq.  (20)  upon  the 
equates  oT  the  ajB.  In  the  limit  of  a  continuous  spectrum,  the  right  hand 
side  can  be  shown  to  approach  aero. 

It  should  be  noted  that  there  are  other  causes  for  tension  variations, 
in  particular  the  creation  of  a  longitudinal  wave  Us  the  transverse  wave 
travels  down  the  wire.  This  is  a  first-order  effect,  however,  and  it  will 
be  taken  into  account  in  the  later  sections  of  this  report. 
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Ill  EXACT  DESCRIPTION  OF  THE  WIRE  MOTION— 
FIRST-ORDER  THEORY 


A.  COORDINATE  SYSTEM 

The  preceding  chapter  haa  given  an  approximate  theory  of  tranaverae 
wave  motion  on  the  trailing  wire,  baaed  on  the  aaaumptiona  that  there  ia 
no  longitudinal  motion  and  that  the  tenaion  in  the  wire  doea  not  vary  with 
time.  A  more  complete  theory, 
including  theae  effects,  will 
now  be  given.  y 

The  position  of  a  point  on 
the  cable  will  be  described  in 
terms  of  the  Cartesian  coordi-  * 

nates  of  its  equilibrium 
position,  and  yQt  and  its 
displacements  from  equilibrium 
in  directions  normal  and 
tangential  to  the  cable  at  that 
point,  n  and  p  respectively 
(See  Fig.  4).  The  ordinary 
cartesian  coordinates  of  the 

point  are  related  as  follows  COORDINATES  OF  A  POINT  ON  THE  WIRE 

to  those  just  described: 

x  -  xQ  +  p  cos  <£0  -  n  sin  <p0  (21a) 

y  c  y0  +  P  **n  +  n  cos  *0  •  (21b) 

B.  EQUATIONS  OF  MOTION 

An  element  of  the  wire  in  vibration  is  acted  upon  as  before  by  the 
forces  shown  in  Fig.  1.  Now,  however,  the  forces  and  angles  may  be 
functions  of  time  and  dT  is  no  longer  identically  zero.  The  component 
of  air  velocity  normal  to  the  element  of  wire  is  given  by 
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(22) 


3n 

r.  .  v, 

where  <f>  is  the  actual  (displaced)  direction  of  the  element.  If 
'bn/'bt  «  V#  sin  <£0,  then  fro*  Eq.  (2)  one  can  write  for  the  aerodynamic 
force 


pV^  /  2  dn  \ 

dF  *  CB  — —  D  ds  ain2  ^  “  *JT  ^ 


Fro*  Fig.  1  it  can  be  seen  that 


dF  +  T  ~  d,  «  pdt  - 

«  3t2 


(23) 


(24) 


where  n  represents  the  cable  mass  per  unit  length.  Now  let 
T  «  (1  +  /3) ,  where  Tp  is  the  average  tension ,  and  fi  is  a  function  of 

position  and  ti»e,  assumed  ssiall  compared  to  unity  for  the  purposes  of 
this  theory.  Let  V 0  *  /T0/m  represent  the  ordinary  wave  velocity 
corresponding  to  the  average  tension.  In  addition,  the  following  new 
dimensionless  quantities  are  defined: 


N 


n 

C 


P 


£ 

C 


r 


v 

c 


The  quantities  C,  77,  &,  and  y  are  defined  as  in  Sec.  II  of  this  report. 
Eq.  (24)  now  becomes 


(i  ' 27  T^ctc  *) 


sin*  4>  +  (1  +  /3)  — *  *  * 
do¬ 


lt 


dr* 


(25) 


For  the  static  esse  Bty/Br* ,  BlV/Br,  and  y9  are  all  equal  to  aero,  and 
Eq.  (25)  gives  the  static  shape  of  the  cable:  sin2  4>0  +  Btf>0/Bcr  *  0,  and 
a  *  cot  <p$t  as  before.  It  is  desirable  to  expand  <f>  about  as  follows: 


<£(ct,t)  *  <fi0(<r)  +  a(crfT) 


(26) 


It  is  possible  to  show,  by  weans  of  Eqs.  (21a),  (21b),  and  (26),  and 
elementary  trigonometry,  that 


a 


—  +  P - 

Btr  Bcr 


(27) 


With  a  «  1, 


sin*  <f>  *  sin*  <f> 0  +  2a  sin  <f>0  cos  4>0  (28a) 

esc  4>  *  esc  4>q  -  a  esc  cot  <t>9  .  (28b) 

Eqs.  (26),  (27),  (28a)  and  (28b)  may  now  be  substituted  into  (25), 
a  and  /3  being  treated  as  small  quantities  whose  products  and  squares  are 
neglected.  Using  the  known  equilibrium  shape  of  the  wire,  and  remembering 
from  Eq.  (22)  that  y  (B/V/Br)  «  sin  <f>,  one  obtains  finally  the  following 
differential  equation  of  motion  normal  to  the  cable: 

a**  2tr  av  i  bp  i  „  b*jv  i  a/v 

- +  -  —  -  - — - fi  »  - +  Try  - — n  —  . 

Bo*  l+tr1^  l  +  a*  ^°r  1  +  a*  Br*  ( 1  +  o* ) 

(29) 

A  second  examination  of  Fig.  1  leads  to  the  following: 


BT  B*p 

— -  ds  ■  jjds  -  (30) 

Bt* 


or 


3£_  B*P 
^  ’  Br* 


(31) 


This  is  the  differential  -equation  for  motion  tangential  to  the  cable. 


C.  INEXTENSIBILITY  OF  THE  CABLE 


Eqs.  (29)  and  (31)  give  two  relations  among  the  three  variable  N,  P , 
and  )9,  The  problem  ia  therefore  not  yet  sufficiently  determined  for  a 
solution.  A  third  relation  may  be  obtained  from  the  elastic  properties 
of  the  cable.  If,  as  before,  the  cable  is  considered  to  be  inextensible , 
one  must  have 


dx2  *  dx2  +  dy2  ■  dx2  +  dyj  .  (32) 


With  the  help  of  Eq.  (21a, b),  and  remembering  that  n  and  p  are  small 
compared  with  xQ  and  y0,  one  easily  obtains 


dp 

n  IT  *  "b7 


or 


N 


-(1  +  cr2) 


dcr 


(33) 


(34) 


Eq.  (34),  together  with  Eqs .  (29)  and  (31),  yields  in  principle  a  complete 
solution  to  the  problem,  provided  the  appropriate  boundary  conditions  are 
specified. 

D.  SIMULTANEOUS  SOLUTION  OF  PARTIAL  DIFFERENTIAL  EQUATIONS 

It  is  possible  to  eliminate  either  P  and  J3  or  N.  and  from  the  system 
of  three  partial  differential  equations.  The  result  in  either  instance  is 
a  fourth-order  linear  partial  differential  equation  whose  coefficients  are 
rather  complicated  functions  of  cr. 

Because  the  coefficients  are  independent  of  time,  it  is  possible  to 
find  product  solutions  in  whi^h  the  time  occurs  only  through  the  factor 
€*v*~i.e.,  let 


N(cr,T )  * 

»(cr)e iVT 

(35a) 

PW,r)  - 

v(a)eiyT 

(35b) 

yS(a,r)  » 

r(cr)  eivT 

(35c) 
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Such  solutions  represent  the  steady-state  behavior  of  the  wire  under  the 
influence  of  a  single  sinusoidal  excitation  of  dimensionless  frequency  v, 
defined  as  follows: 


(36> 


where  0)  is  the  angular  frequency  of  the  input. 

Formal  substitution  of  Eq.  (35a)  or  Eq.  (35b)  into  one  of  the  above- 
mentioned  fourth-order  equations  converts  it  into  a  linear,  fourth-order, 
ordinary  differential  equation  with  complex,  non-constant  coefficients. 

In  principle  it  is  then  possible  to  obtain  the  functions  w(cr),  v(cr)t 
r(a)  in  a  particular  case,  provided  four  suitable  boundary  conditions  are 
supplied.  In  practice,  however,  the  nature  of  the  differential  equations 
is  such  that  it  is  nearly  impossible  to  obtain  any  sort  of  general 
information  from  them.  Tin  fact,  any  attempt  to  find  approximate  forms  of 
a  literal  solution  becomes  rapidly  so  involved  as  to  obscure  entirely  the 
bas  ic  phenomena  which  it  is  desired  to  represent. 


E.  THE  UNIFORM-TENSION  CASE 

In  many  cases  of  physical  interest,  the  total  mass  of  the  wire  is 
very  much  smaller  than  that  of  the  end-weight.  Under  these  circumstances , 
the  fluctuation  in  tension  produced  at  the  lower  end  of  the  wire  by 
acceleration  of  the  end-weight  is  much  larger  than  any  variation  with 
distance  along  the  wire  produced  by  longitudinal  acceleration  of  the  wire 
itself.  This  is  especially  true  for  long  (large  Z)  wires  where  the  inputs 
are  largely  transverse,  as  then  any  longitudinal  motion  is  confined  almost 
entirely  to  the  lower  part  of  the  wire,  between  a  «  0  and  <r  -  2. 

If  the  conditions  described  above  are  satisfied,  the  tension  can  be 
considered  to  be  independent  of  or,  although  the  time  dependence  remains. 

In  other  words,  the  tension  variation  amplitude  r(cr)  becomes  a  (complex) 
constant,  whose  value  is  determined  by  the  boundary  conditions  at  cr  «  0. 
Eqs.  (34),  (35a),  and  (35c)  may  now  be  substituted  into  Eq.  (29),  giving 

w"  +  — -  v*  +  f- - - -  +  v2  -  2iyv  - - — rr\  w  *  - ~ -  r 

1  +  cr1  lil  +  a2)  (1  +  c2)  J  1  +  a2 

(37) 
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where  primes  denote  differentiation  with  respect  to  a. 

A  further  simplification  may  be  effected  through  the  following 
substitution : 

w(cr)  «  (1  +  cr2 )  u(cr)  .  (38) 

The  result  is 

u"  +  v*  /(cr)  u  *  (i  +  a2)  ^  r  ,  (39) 

where 

f(cr)  =  1  -  2i(y/v){\  +  cr2)  ^  .  (40) 

The  general  solution  of  Eq.  (39)  will  consist  of  a  homogeneous  solution 
containing  two  arbitrary  constants,  plus  a  particular  integral  which  will 
be  proportional  to  r. 

If  /(cr)  is  sufficiently  slowly-varying,  a  very  good  approximation  to 
the  homogeneous  solutions  is  given  by 

u  «  const  x  exp[±i  v  f  Sf(cr)  da]  .  (41) 

For  the  range  of  values  of  y  and  v  represented  by  the  shaded  portion  of 
Fig.  5  (Regions  I  and  II),  the  approximation  of  Eq.  (41)  holds,  and  it 
is  possible  to  obtain  analytic  solutions  of  the  following  form: 

Uj(a)  «  (cr  +  /l  +  a2)7  exp{i  v  [cr  +  (1  /2)(y/v)2  tan-*1  cr] } 

( D-14a) 

Uj  (cr)  •*  (cr  +  /l  +  cr*  )  y  exp{-i  v  [cr  +  (1/2 )(j/v)2  tan"1  cr]^ 

(D-14b) 

Details  of  the  derivation  of  (D-14a,b)  and  the  determination  of  the 
boundaries  in  Fig.  5  may  be  found  in  Appendix  D. 

The  arc-tangent  term  represents  a  deviation  of  the  wave  velocity 
from  its  normal  value  U0  which  is  not  sufficient  to  alter  the  amplification 
in  this  order  of  approximation.  This  term  can  be  ignored  in  Region  I  of 
Fig.  5,  but  even  in  Region  II,  where  strictly  speaking  it  cannot,  only 
the  phase  of  the  wave  at  points  other  than  the  origin  is  affected. 
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FIG.  5 

DOMAINS  OF  VALIDITY  OF  SOLUTIONS 

For  the  sake  of  mathematical  simplicity,  the  remainder  of  this 
investigation  will  be  limited  to  the  case  where  this  term  is  neglected , 
with  the  understanding  that  it  could  be  included  if  necessary.  Using 
Eqs.  (35a)  and  (38),  one  then  has 

Nx(ct,t)  -  (l+cra)^(cr  +  /l  Va*)y  elv(r4tr)  (42a) 

-  (1  +  a*)"54  (a  +  A  +  aJ)‘r  e**™*  (42b) 

which  are  easily  seen  to  be  identical  with  the  solutions  for  the  transverse 
wave  obtained  by  the  simple  approach  of  Sec.  II  of  this  report. 
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It  will  be  remembered  that  the  complete  solution  of  Eq.  (39)  consists 
of  the  homogeneous  solutions,  Eqs.  (D-14a,b)  plus  s  particular  integral. 
The  latter  may  be  represented  as  r  u$(a),  where  u3(<r)  i*  *  particular 
integral  of  the  following  equation: 

u'  +  v2f(cr)  a  *  (1  +  or1)  .  (43) 

There  are  several  well-known  and  straightforward  methods  for  finding 
an  integral  of  such  an  equation.  In  this  case,  however,  all  require  one 
or  more  numerical  integrations  for  each  value  of  v  which  it  is  desired  to 
consider.  Also,  since  there  is  an  infinite  number  of  particular  integrals 
differing  by  arbitrary  multiples  of  the  homogeneous  solutions,  there  is 
no  assurance  that  one  of  these  methods  will  yield  an  integral  which  will 
aid  in  the  physical  interpretation  of  the  problem. 

If  the  range  of  cr  which  is  of  interest  in  ^  particular  case  is  known, 
a  suitable  approximate  expression  for  (1  +  cr2 )  may  be  found  in  a  form 
which  can  be  handled  analytically.  An  expression  in  the  form  of  a  sum 
of  exponentials  is  especially  useful,  as  it  leads  to  an  integral  which 
becomes  very  small  for  values  of  cr  greater  than  a  few  times  unity  — i.e., 
for  the  less  sharply  curved  portions  of  the  wire.  In  addition,  for  many 
cases  the  solution  is  very  nearly  pure  real  and  positive,  corresponding 
to  a  sort  of  heaving  motion  of  the  wire  as  a  whole.  (See  Sec.  IV  for  a 
specific  application  of  this  approach.) 

Thus,  one  is  led  to  a  representation  of  the  wire  motion  in  terms  of 
a  downward-traveling  wave  uJf  a  reflected,  or  upward-traveling,  wave  Uj , 
and  a  heaving  motion  u3  which  is  most  important  in  the  lower,  more 
sharply  curved,  portions  of  the  wire.  The  general  motion  of  the  wire 
can  be  written  as 

N(cr,r)  *  (1  +  o3)  ^  [fcjUjta)  +  b2uJ(cr)  +  r  Uj(cr)]  eivr  (44) 

where  61  and  are  arbitrary  constants  yet  to  be  determined  and  r  is  the 
complex  amplitude  of  the  fractional  tension  fluctuation  in  the  wire. 

The  general  longitudinal  motion  of  the  wire  can  be  obtained  from 
Eq.  (44)  by  means  of  Eq.  (34).  The  result  is 
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P(cr,r) 


(45) 


»  [t>(0)  -  bjVjtcr)  -  5jVj (<r)  -  r  i/,(a)]  eivr 

where 


vn(ar) 


(1  +  cr 


'*)  ^  u.(a!)d<r' 


(46) 


F,  APPLICATION  OF  BOUNDARY  CONDITIONS 

It  will  be  remembered  that  the  exact  solution  of  the  problem  revolves 
about  a  fourth-order  linear  differential  equation.  Since  such  an 
equation  requires  the  specification  of  four  boundary  conditions,  any 
approximate  solution,  to  be  valid,  must  contain  the  same  number  of 
conditions. 

If  the  case  of  interest  is  that  in  which  the  input  is  at  the  upper 
end  of  the  trailing  wire,  then  the  first  two  conditions  will  be  the 
specification  of  JV(2,r)  and  P(£,r),  the  transverse  and  longitudinal 
inputs  at  the  upper  end  of  the  cable.  The  remaining  two  conditions  will 
depend  upon  the  manner  of  termination  of  the  cable  at  its  lower  end. 

One  will  be  a  relation  connecting  y8(r) ,  the  tension  fluctuation,  with 
P(0,r),  the  longitudinal  response  of  the  terminating  body.  The  other 
will  describe  the  transverse  response  of  the  same  body  by  means  of  a 
relation  between  cW/3ct(0,t)  and  iV(0,r);  Now  N(2tr)  »  w(Z)elVr,  and  an 
application  of  Eq.  (44)  at  <r  -  £  yields 

(1  +  I*)X  .(Z)  -  +  6jUj(Z)  +  r  u,(Z)  .  (47) 

Proceeding  in  the  same  manner  with  Eq.  (45),  one  obtains 

v(Z)  *  v(0)  -  blvl{l)  -  b2v2(l)  -  r  Vj(2)  .  (48) 

The  varying  tension  in  the  cable  will,  in  general,  produce  a  response 
of  the  terminating  weight,  which  will  depend  upon  the  method  of  attach¬ 
ment,  damping  forces  present,  etc.  This  response  can  always  be  expressed 
by  means  of  a  complex  coupling  constant  £ p ,  as  follows: 

v(0)  -  .  (49) 
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Similarly,  there  will  be  a  transverse  motion  of  the  weight  in  response 
to  the  transverse  component  of  the  tension  force  applied  to  it: 

*<0)  «  £,*'<<»  .  (50) 

This  is  equivalent  to  the  following: 

+  bt  +  r  Uj (0)  *  (y  +  iv)£Mbl  -  (y  +  iv)Cjr^2  +  Cg  r  "j(0) 

(51) 

If  Eq.  (49)  is  substituted  into  Eq.  (48),  and  Eqs.  (47)  and  (51) 
are  rewritten,  one  obtains  the  following  system  of  algebraic  equations: 

Ui (X>61  +  u2(Z)62  +  u3(2)r  «  (1  +  T2)^  w(2) 

t fxa)bx  +  v2(2)62  +  iv%(2)  -  t,]r  -  — w(I)  1(52) 

[1  -  (7  +  iv)£j9)bl  +  [1  +  (7+  iv)C^]62  +  [u3(0)  ^£tt»,s(0)]r  -  0  . 

Eqs.  (52)  can  readily  be  solved  for  6lf  6a,  and  r,  if  the  four  boundary 
conditions  are  expressed  in  terms  of  w(Z),  v’(Z),  £p ,  and 

In  addition  to  the  vibrational  inputs  expected  at  the  point  of 
attachment  to  the  airframe,  there  will  be,  depending  upon  speed  and 
details  of  the  mechanical  configuration,  other  inputs  at  points  farther 
down  the  wire  resulting  from  the  action  of  oblique  shock  waves  originating 
on  the  airframe.  Because  of  their  lower  point  of  application,  the 
amplification  will  be  less,  but  the  net  effect  is  difficult  to  assess 
theoretically  in  the  absence  of  experimental  data. 

G.  LOW  FREQUENCY  BEHAVIOR 

It  is  of  considerable  interest  to  examine  the  behavior  of  low- 
frequency  waves  on  the  wire  under  conditions  for  which  the  approximations 
of  the  preceding  theory  do  not  hold  (Region  III  of  Fig.  5).  It  will  be 
remembered  that  the  approximate  theory  of  Sec.  II  gives  correct  results 
for  the  transverse  motion  for  the  case  where  the  end  is  longitudinally 
free— that  is,  where  the  longitudinal  motion  is  ignored.  More  importantly, 
it  gives  the  correct  value  for  the  amplification  undergone  by  the 


2( 


downward- t rare ling  wave,  which  is  the  primary  factor  governing  the 
amplitude  of  transverse  notion  near  the  lower  end  of  the  wire  and  the 
severity  of  the  tension  fluctuations  if  the  end  is  restrained. 

If  the  propagation  of  a  wave  in  a  straight  piece  of  wire  is 
considered,  it  is  easily  found  that,  for  sufficiently  high  frequencies, 
the  wave  velocity  is  independent  of  frequency  and  of  the  angle  nade  with 
the  airstream,  having  a  value  of  However,  for  lower  frequencies, 

the  velocity  depends  on  the  angle  and  in  general  increases  with  decreasing 
frequency.  The  ratio  of  this  velocity  to  the  high  frequency  value 
(c/.  W.  H.  Phillips5)  is  given  by 

—  ■  W2  |(v2  sin2  4>  cos2  4>)  +  [^v2  -  sin2  4>  cos2  <£)2  +  4yjv2  sin  <^J  ^ 

(53) 

This  expression  nay  be  used  to  compute,  for  any  given  v,  a  corrected 
value  of  y.  Then,  ignoring  the  effect  of  wire  curvature,  Eq ,  (13)  leads  to 


[cos  <(>  -  y(4>tv)]  CSC  4>  d4> J  (54) 

This  amplification  is  shown  as  a  function  of  v  in  Fig.  6  for  a  wire  of 
£  *  10  and  the  flight  parameters  of  Sec.  IV  of  this  report  —  i.e.,  for 

*  0.23.  The  results  can  be  interpreted  to  indicate,  at  least  qualita¬ 
tively,  that  the  amplification  decreases  rapidly  for  very  low  frequencies. 

However,  it  is  also  found  that  the  attenuation  of  the  reflected  wave 
computed  in  this  way  decreases  to  unity  at  these  low  frequencies.  As  a 
result,  there  is  little  net  aerodynamic  damping  in  the  path  from  the  upper 
end  of  the  cable  to  the  lower  and  back  again,  and  the  resonances  which 
occur  for  certain  frequencies  can  therefore  be  very  extreme  unless  there 
are  other  sources  of  mechanical  damping. 
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FIG.  6 

INCIDENT  WAVE  AMPLIFICATION  AT  LOW  FREQUENCIES 


H.  NON*  AERODYNAMIC  DAMPING 

The  analysis  to  this  point  has  entirely  neglected  damping  froa  sources 
other  than  the  airstrean  itself.  Such  damping  nay  arise  through  the  action 
of  viscous  forces  within  the  cable,  or  through  the  transfer  of  energy  to 
another,  coupled,  node  of  oscillation. 

The  first  type  of  danping  can  be  described  in  terns  of  a  danping 
constant  K  characteristic  of  the  particular  cable  used.  The  effect  will 
be  to  decrease  the  anpli fication  of  the  incident  wave  by  a  factor  exp(~XC£), 
and  to  increase  the  danping  of  the  reflected  wave  by  the  sane  factor. 
Qualitatively,  the  anplitudes  near  the  lower  end  will  be  reduced  for  any 
set  of  conditions,  while  the  resonant  fluctuations  of  these  anplitudes  with 
cable  length  will  becone  less  pronounced  because  of  the  decrease  in  re¬ 
flected  wave  anplitude  at  the  input  point  (the  upper  end)  relative  to  the 
incident  wave. 
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The  quantitative  effect  nay  perhaps  Most  easily  be  seen  by  rewriting 
Eqs.  (19)  and  (42a, b)  to  include  this  daaping: 

«j  -  —  nj  /l  +  Z*  c.ch  [y  .inh"1  I  +  KfS.  +  ivi]  (19) 

2 

/V,(fr,T)  -«  (1  +  (cr  +  /I  +  c r*)r  e*c<r-e  *v<T+<r>  (42*) 

*t(e,r)  -  (1  +  <r*)”X  (o-  +  /I  +  <r*)'y  €'*<*  eiv(r~a>  .  (42b) 

The  nost  inportant  case  of  the  second  type  is  the  coupling  which  is 
certain  to  occur  between  the  oscillations  in  the  vertical  plane  and  lateral 
oscillations  perpendicular  to  that  plane.  As  pointed  out  in  Sec.  II-B, 
the  latter  kind  of  notion  experiences  only  slight  aerodynamic  amplification 
or  damping. 

Lateral  vibrations  will  arise  in  two  ways.  First,  there  will  be 
direct  lateral  inputs  at  the  upper  end  of  the  wire  and  to  a  leaser  extent 
along  the  wire.  Second,  there  will  be  the  coupled  energy  from  the  in¬ 
plane  vibrations,  as  mentioned  above.  The  total  attenuation  in  a  round 
trip  to  the  lower  end  of  the  cable  and  back  to  the  input  point  is  actually 
less  for  these  waves  than  the  net  attenuation  for  the  in-plane  waves. 
However,  because  there  is  no  amplification  of  the  incident  waves,  their 
contribution  to  the  amplitude  at  the  lower  end  will  be  small.  Without 
detailed  knowledge  of  the  coupling  mechanism  it  is  most  difficult  to 
analyse  this  effect  in  a  quantitative  way,  but  qualitatively  it  can  be 
seen  that  the  total  amplitudes  on  the  lower  parts  of  the  wire  will  be 
reduced. 

Both  of  the  damping  mechanisms  described  above  will  be  more  effective 
for  a  stranded  cable,  in  the  one  case  because  of  the  obvious  increase  in 
internal  friction,  in  the  other  because  the  helical  form  of  individual 
wire  strands  tends  to  rotate  the  plane  of  vibration  of  the  wire,  converting 
an  in-plane  vibration  to  a  lateral  vibration. 
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IV  APPLICATION  OF  THE  THEORY  TO  A  PARTICULAR  PROBLEM 


A.  FLIGHT  PARAMETERS  AND  CONSTANTS  OF  AN  EXPERIMENTAL 
TRAILING- WIRE  ANTENNA 

Computation*  have  been  made  for  an  experimental  trailing- wire  antenna 
whose  mechanical  properties  have  been  tested  in  flight  at  about  35t000  feet 
altitude  and  speeds  up  to  Mach  1,24.  Results  of  these  highly  successful 
flight  tests  will  be  presented  as  a  part  of  the  final  report  on  another 
SRI  project.* 

The  cable,  1/16  inch  in  diameter,  was  tested  at  lengths  up  to  150  feet. 
It  is  terminated  in  a  19-pound  f in- stabili xed  weight  of  high  fineness  ratio, 
attached  at  its  center  of  gravity  to  the  cable.  Two  systems  of  springing 
were  investigated,  one  an  essentially  rigid  coupling,  the  other  a  pre¬ 
stressed  spring  permitting  approximately  a  i  10-percent  variation  in  cable 
tension  with  s  1  1-inch  travel  in  the  vertical  direction.  The  drag  on  the 
edd-weight  in  flight  is  quite  small,  so  that  the  cable  is  essentially  ver¬ 
tical  at  its  lower  end,  and  the  static  value  of  tension  is  very  nearly 
equal  to  the  weight,  or  19  pounds. 

The  accepted  value  of  the  drag  coefficient  for  a  circular  cylinder  of 
infinite  length,  at  moderate  subsonic  Mach  numbers  and  subcritical  Reynolds 
numbers,  is  1.20.  As  the  Mach  number  is  increased  above  0.3,  the  coef¬ 
ficient  riaes  with  some  irregular  variations  to  a  maximum  of  2.13  just  below 
sonic  speed.  Beyond  this  maximum,  a  decline  occurs  toward  an  asymptotic 
value  of  1.33. 12  Clearly,  the  value  of  drag  coefficient  for  a  trailing-wire 
antenna  will  be  a  function  of  the  speed  of  the  aircraft,  and  also  will  vary 
with  position  along  the  cable.  However,  the  physically  interesting  case  of 
high  wave  amplification  always  implies  large  values  of  dimensionless  cable 
length  Z,  and  therefore  very  small  angles  of  attack  over  most  of  tbe  length 
of  the  cable  (see  Fig.  2).  As  a  result,  even  when  the  aircraft  speed  is 
highly  supersonic  tth'e  transverse,  airflow  across  most  of  the  wire  is 
subsonic.  It  is  possible  to  compute  bothtthe  cable  shape  and  the  total 
amplification  exactly  in  a  given  case  by  means  of  numerical  integration. 


*  SRI  Project  1404,  " InTeatifa tion  if  ia  Excited  Airfraee  ee  ai  Aateaaa.” 
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However,  it  hn  been  found  that  a  conatant  value  of  C D  ■  1.25  given  the 
correct  value  of  amplification  within  about  5  percent  over  a  very  wide 
range  of  flight  conditions,  and  it  haa  been  adopted  for  the  sake  of  com* 
putational  simplicity.  On  this  basis,  and  with  the  additional  information 
that  \i%  the  cable  mass  per  unit  length,  is  2,34  *  10=4  slug/ft,  one  finds 
for  Mach  1.24  at  35,000  feet  altitude  and  a  typical  cable  length  of  100  feet* 

7  •  0.23 
C  »  5.25  feet 

I  -  19.0  . 

The  cable  shape  for  these  conditions,  on  the  other  hand,  corresponds  more 
nearly  to  C D  ■  1.44. 

The  weight  of  the  wire  at  this  length  is  only  1.13  pounds,  and  there¬ 
fore  the  uni  form- tension  approximation  of  Sec.  III-E  is  justified.  Ref¬ 
erence  to  Fig.  5  shows  that  the  system  parameters  lie  in  Region  I,  where 
the  simplest  form  of  solution  holds  if  v  is  greater  than  about  0.88; 
corresponding  to  an  actual  frequency  of  7.6  cycles.  The  aerodynamic 
amplification  of  the  downwind  wave  is  8.2. 

The  system  described  above  is  designed  to  simulate  the  aerodynamic 
amplification  of  a  hypothetical  higher-speed  system  identical  to  it  in  all 
respects  except  for  the  use  of  a  50-pound  end-weight.  On  the  basis  of  a 
100- foot  wire,  the  following  are  typical  speed- altitude  combinations  for 
thelhypothstical system  which  give  the  same  theoretical  amplification  as  the 
experimental  system  does  in  the  Mach  1.24  flight  test: 

Mach  2.0  -  35,000  feet 
Mach  2.5  -  48,500  feet 
Mach  3.0  -  59,000  feet. 

It  should  be  noted  that,  although  the  wave  amplification  is  the  same  in 
all  the  above  cases,  the  cable  shape  differs  from  one  case  to  another  be¬ 
cause  each  represents  a  different  y,Z  combination. 
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B.  COMPUTATION  OF  THE  uH(a)  AND  v  n(<?) 

The  magnitudes  of  the  two  homogeneous  solutions  u^o-)  and  u^o*) 
(downwind  and  upwind  wares,  respectively)  depend  only  upon  y,  and  hare 
been  shown  in  Fig.  7  for  y  ■  0.23.  The  phases,  of  course,  depend  upon 
position  through  the  factors  etiv<r . 

The  particular  integral  u3(<x)  is  a  solution  of  Eq.  (43),  which  is 
here  repeated  for  convenience: 


u"  +  v%  f(cr )  u  *  (1  +  a2)  ^ 


y  '  0.23 


FIG.  7 

MAGNITUDES  OF  TRAVELING  WAVE  SOLUTIONS 


The  solution,  which  is  in  general  complex,  will  be  written  in  the  following 
form: 

us(<t)  «  x(cr)  +  i  2r/v  y(cr)  (55) 

where  x(<r)  and  y(cr)  are  real.  If  Eq.  (55)  ia  substituted  into  Eq.  (43), 
and  real  and  imaginary  parts  of  the  resulting  equation  are  written  sepa¬ 
rately,  one  obtains  the  following: 

x"  +  v2  [pc  +  4y2  (1  +  or2)  ^  y]  «  (1  +  cr2)  ^  ,  (56a) 

y"  +  v2  y  *  (1  +  a2)  ^  x  .  (56b) 

This  system  of  equations  can  be  solved  quite  readily  for  x  and  y  if 

)4y2  ( 1  +  or2)  ^  y  |«|x|  •  (57) 

A  procedure  for  finding  a  particular  integral  of  Eqs.  (56a, b)  under  the 
assumption  that  Eq.  (57)  holds  is  given  in  Appendix  E.  Resulting  values 
of  u3(cr)  are  shown  in  Figs.  8(a),  (b)  for  the  experimental  system  at 
several  frequencies. 

Detailed  examination  of  a  number  of  cases  shows  that  at  worst 
4  (v  ^  1,  cr  •  0)  the  left  side  of  Eq.  (57)  is  about  13  percent  of  the  right 
side.  Although  this  would  seem  to  make  the  approximation  rather  marginal, 
both  |y  |x  |  and  (1  +  cr2)  decrease  rapidly  with  cr,  and  \y/x\  also  decreases 
with  v.  Thus  the  neglected  term  is  nearly  always  much  smaller  than 
13  percent  of  the  leading  term. 

The  longitudinal  wave  functions  vn(cr)  can  readily  be  obtained  from 
the  un(cr)  through  numerical  evaluation  of  Eq.  (46).  The  integral  converges 
very  rapidly  to  a  constant  value  v^I00)  for  cr  £  5.  These  asymptotic  values 
are  shown  in  Table  I. 


TABLE  I 


V 

1 

2 

3 

4 

5 

0.851  /  41.4° 

0.586  /  63. 1° 

0.414  /  76.2° 

0.305  /  83.9* 

0.234  /  88. 1® 

v2(v>) 

0.689  /-34.80 

0.508  /— 55. 4° 

0.378  /-68.60 

0.288  /-76. 9® 

0.227  /-82.0® 

*,<«> 

0.652  /  12.9° 

0.1915  /  9.1° 

0.0874  /  6.6° 

0.0493  /  5.1® 

0.0315  /  4.2® 
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FIG.  8(b) 
PHASE  OF  u3 


C.  RIGID  TERMINATION 

If  the  lower  end  of  the  cable  is  attached  directly  to  the  center  of 
mass  of  a  heavy,  stable  end-weight,  practically  no  motion  of  that  point 
will  take  place.  Mathematically,  this  condition  is  expressed  as  follows 
in  terms  of  the  two  mechanical  coupling  coefficients: 

m  0  ,  iH  *  0  (58a, b) 

Let  us  suppose  that  the  vibrational  input  at  the  upper  end  of  the  wire  is 
purely  transverse  and  of  unit  amplitude — i.e., 

v(X)  -  0  .  •»(!)  -  1  (59a, b) 

Further  let  it  be  assumed  that 


2  Z  5 

36 


Then  the  amplitudes  of  the  incident  and  reflected  wares,  os  well  as 
the  resulting  tension  fluctuation,  can  be  computed  with  the  aid  of 
Eqa.  (52),  which  may  now  be  written  as  follows; 


Ka) 

«,(£) 

V  J  (00) 

v®) 

\  1 

1 

(61) 


Solutions  for  6 2 ,  and  r  may  be  obtained  in  the  following  form; 


B  A  +  2* 

bx  ■  - — -  (62a) 

B.  A  +  Is 

f>2  "  - (62b) 

Bj  A  +  2J 

r  =  (62c) 

D 

where 

*  Vj(oo)  u3(0)  "  t^f00)  (63a) 

B2  *  -[v^co)  u3 ( 0)  -  *3(“>)]  (63b) 

B3  -  vt(oo)  -  t/2(®)  (63c) 

D  -  B^^S)  +  B,u2(2)  +  B3u3(Z)  (63d) 


The  absolute  magnitudes  of  the  solutions,  Eqs.  (62a, b,c,),  pass  through 
maxima  and  minima  as  2  is  raried,  corresponding  to  resonances  in  the 
standing  ware  on  the  wire.  The  principal  factor  controlling  the  location 
and  height  of  these  maxima  is  the  quantity  |B|. 

The  complex  factors  Bn  are  constant  for  these  relatirely  large  ralues 
of  2.  In  describing  the  total  transrerse  motion  of  the  wire  according  to 
Eq.  (44),  the  term  6j  u^cr)  can  be  represented  by  a  rector  which  increases 
slowly  in  length  with  cr  while  rotating  counterclockwise  at  a  rate  of 
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v  radians  per  unit  of  cr;  the  vector  representation  of  63  Uj(a)  decreases 
slowly  with  cr  while  rotating  clockwise  at  the  same  rate;  and  r  u5(<r) 
decreases  quite  rapidly  with  cr,  maintaining  very  nearly  a  constant  angular 
position. 

A  numerical  evaluation  of  the  various  quantities  in  Eqs.  (63)  shows 
that  the  third  term  in  D  is  always  considerably  smaller  than  either  of 
the  other  terms.  Thus,  the  magnitude  of  the  resultant  D  will  be  a  minimum— 
i.  e.,  the  response  will  be  a  maximum— when  the  first  two  terms  are  very 
nearly  antiparallel.  On  the  other  hand,  minimum  responses  (maximum  |/)|) 
will  occur  when  the  two  dominant  terms  are  nearly  in  phase. 

The  angular  position  of  the  third  term  is  such  that,  in  the  minimum 

response  condition,  it  is  nearly  parallel  or  antiparallel  to  the  resultant 

of  the  first  two  terms.  As  a  result,  there  are  two  classes  of  minima, 
which  alternate  as  X  varies,  one  lying  on  a  higher  envelope  than  the  other. 
In  the  maximum  response  condition,  however,  the  third  term  has  practically 
no  effect  on  the  magnitude  of  D. 

Of  the  greatest  physical  interest  is  the  “worst”  case  (maximum 
response),  which  may  be  obtained  by  ignoring  the  third  term  in  D  entirely 

and  treating  the  first  two  terms  as  antiparallel.  Values  of  blt  b2,  and 

r  have  been  computed  in  this  way  for  various  values  of  X  and  v ,  on  the 
assumption  that  each  input  point  is  a  node  in  the  standing  wave,  regardless 
of  the  actual  phase  relationships  at  the  point  in  question.  The  results 
have  been  shown  in  Tables  II  and  III.  In  addition,  Fig.  9(a),  (b),  (c), 

(d)  shows  the  way  in  which  the  three  waves  combine  at  the  origin. 
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D.  TERMINATION  WITH  ZERO  SPRING  CONSTANT 

If  one  postulates  a  cable  termination  consisting  of  a  pre- stressed 
spring  giving  the  proper  value  of  equilibrium  tension,  but  permitting 
longitudinal  motion  with  a  dynamic  spring  constant  of  zero,  the  boundary 
conditions  are  expressed  by  «  0,  r  •  0,  t>(2)  ■  0,  and  v(Z)  ■  1/0°. 
Although  the  tension  variation  r  is  zero,  the  longitudinal  motion  of  the 
lower  end  of  the  cable,  WO)  [  ■  Cpr],  is  finite.  For  this  case,  Eqs.  (78) 
may  be  written 


/u,(Z)  uj(I)  °Wbi\ 

«v®)  *a<®)  -1  |  b2 

V  1  1  0 J  \v(0 )/ 


/a 


TT5\ 


(64) 


) 


The  first  and  third  components  of  Eq.  (64)  can  be  solved  independently  of 
the  second,  since  Eq.  (64)  can  also  be  written  as  follows: 


/ul<z) 

u2(I)  \ 

[M 

fATJi) 

1 

.  1 

Uj 

{  0  / 

wo) 


6 j  Vj(®)  +  b2  v2(cc>) 


(65a) 


(65b) 


The  solution  of  Eqs.  (65a, b)  gives 


A  +  zJ 


WO) 


A  +  zJ 


Uj(Z)  -  u2(Z) 

A"  +  i* 

"  ttj(Z)  -  u2(X) 

«'1(®)  -  V2(C0) 


ui(Z)  ~  u2(l) 


(66a) 


(66b) 


(66c) 


It  is  seen  that  when  the  lower  end  of  the  wire  is  free  to  move  longi¬ 
tudinally,  there  results  a  simpler  type  of  transverse  motion  than  if  it 
is  constrained,  involving  only  u2  and  u2.  Maximum  and  minimum  responses 
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in  this  case  correspond  exactly  to  nodes  and  antinodea  in  the  standing 
wave  pattern  formed  by  the  incident  and  reflected  wares.  Numerical  results 
again  for  the  "worst”  case,  hare  been  obtained  and  are  shown  for  comparison 
in  Tables  II  and  III.  The  two  wares  are  always  equal  in  magnitude  and 
180°  out  of  phase  at  the  origin. 

A  clearer  idea  of  what  happens  on  the  wire  can  be  had  by  plotting 
the  transrerse  amplitude  as  a  function  of  position  along  the  wire.  The 
result,  for  2  ■  10,  v  *  1,  and  for  both  the  fixed  and  free  cases,  is  shown 
in  Fig.  10.  The  amplitude,  for  the  free-end  case,  oscillates  between 
smooth  upper  and  lower  enrelopes,  given  by 

*  6. 77  /l  +  a*  ( |  |  +  |u2(ct)|)  (67a) 

l*’(°r)l.i.  *  6.  77  /l  +  cr2  (  |  Uj  (cr)  |  -  |  u  2  ( cr )  | )  .  (67b) 


FIG.  10 

VIBRATIONAL  AMPLITUDE  AS  A  FUNCTION  OF  POSITION  ON  THE  WIRE  FOR  v  =  l 
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At  higher  frequencies,  the  envelope  regains  the  same,  while  the  fluctuations 
of  amplitude  with  position  become  increasingly  rapid.  The  first  maximum 
is  always  the  highest,  but  can  never  exceed  13.54  tines  the  input  amplitude. 

With  the  rigid  termination,  the  three  solutions  combine  in  a  com  pli¬ 
cated  way  in  the  lower  part  of  the  wire,  yielding  an  amplitude  pattern 
which  does  not  necessarily  stay  within  the  same  two  envelopes.  The  exact 
nature  of  this  pattern  is  dependent  upon  frequency,  and  it  is  possible  for 
the  second  maximum  to  be  as  high  as  about  20  times  the  input  amplitude  at 
some  frequencies.  At  points  farther  up  along  the  cable,  the  standing  wave 
has  much  the  same  appearance  as  in  the  free-end  case,  except  for  a  phase 
shift  of  the  maxima  and  minima  relative  to  their  former  positions,  and  a 
somewhat  greater  fluctuation  in  amplitude  due  to  the  presence  of  a  larger 
proportion  of  than  before.  The  effect  of  decreases  markedly  at 
higher  frequencies,  and  the  whole  pattern  approaches  that  which  is  obtained 
for  the  free  case. 

It  should  be  noted  that,  although  the  transverse  amplitude  patterns 
in  the  two  extreme  cases  are  the  same  at  high  frequencies,  the  longitudinal 
motion  and  tension  fluctuation  are  quite  different.  In  the  free  case  the 
tension,  of  course,  remains  constant,  while  there  is  considerable  longi¬ 
tudinal  motion  of  the  lower  end  of  the  cable,  decreasing  at  the  higher 
frequencies;  in  the  fixed  case  there  is  no  motion  of  the  lower  end,  but 
the  tension  variation  becomes  progressively  more  severe  with  increasing 
frequency. 

E.  OTHER  TERMINATIONS 

It  is  possible  to  consider  many  other  theoretical  terminations,  in¬ 
cluding  viscous  types  capable  of  absorbing  energy.  Since  two  constants, 
in  and  are  needed  to  describe  the  termination  completely,  it  is 

possible,  at  a  given  frequency  and  for  a  particular  set  of  flight  parameters, 
to  satisfy  two  independent  conditions.  It  is  not  difficult  to  show  that 


>KP 


ij/2 

K  HP 


«»>  +  ia 


HP 


(68) 


In  Eq.  (68),  G  represents  the  acceleration  of  gravity  in  consistent 
dimensionless  units,  and  a ^  is  a  dimensionless  coefficient  of  viscous 


/ 
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friction  in  the  M  or  P  direction.  These  quantities  are  defined  aathe- 
natically  as  follows; 


G 


a 


»n9C 


mV 


(69) 


(70) 


where  is  the  coefficient  of  viscous  friction  in  conventional  units — 

e.g. ,  lb-sec/ft.  The  frequency  ^  is  the  (diaensionless)  natural  fre¬ 
quency  of  the  weight  and  its  suspension,  and  m  is  the  Mass  of  the  weight 

For  the  actual  experimental  systen  using  a  pre- stressed  spring,  the 
flight-test  conditions  and  the  constants  of  the  linkage  are  such  that 

G  »  2.07  x  10’3  ,  ® 

vp  £  .0. 13  ,  ap  ^  0 

These  values  lead  to  the  following; 

*  0  ,  ip  ^  0.123^° 

Numerical  results  in  the  coaputation  of  w(a)  are  shown  for  comparison  in 
Tables  II  and  III.  It  is  seen  that  this  termination  is  very  auch  like  the 
fixed  case  at  low  frequencies,  and  very  auch  like  the  free  case  at  high 
frequencies.  It  is  possible  however,  at  certain  frequencies,  to  have 
either  longitudinal  aaplitudes  greater  than  in  the  free  case  or  tension 
fluctuations  greater  than  in  the  fixed  case. 

The  longitudinal  aaplitudes  have  been  shown  as  functions  of  frequency 
in  Fig-  11  for  a  wire  of  ^  *  10  and  fixed  input  amplitude  of  10~3  in 
diaensionless  units,  corresponding  to  0.063  in.  in  this  case.  For  soae 
purposes  it  is  of  greater  interest  to  know  the  longitudinal  response  for 
a  fixed  acceleration.  This  has  been  shown  in  Fig.  12  for  a  diaensionless 
acceleration  amplitude  of  10”*,  corresponding  to  4.82  g.  The  aaplitudes 
of  the  tension  fluctuation,  in  percent,  have  been  shown  for  the  saae  con¬ 
ditions  in  Figs.  13  and  14.  It  is  worthy  of  note  that,  for  aaplitudes  of 
such  size  that  the  cable  never  goes  slack,  the  notion  of  the  lower  end  of 
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FIG,  11 

LONGITUDINAL  AMPLITUDE  AT  LOWER  END  OF  CABLE  WITH  FIXED  INPUT  AMPLITUDE 

the  cable  due  to  the  second-order  effect  described  in  Sec.  II-G  can  be 
neglected  in  comparison  with  the  first-order  effect. 

It  would  be  at  least  theoretically  possible  in  this  same  system  under 
a  particular  set  of  flight  conditions,  to  eliminate  both  the  return  ware 
and  the  “tension  wave"  u}.  First,  the  tension  wave  can  be  eliminated 
by  making  *  00  as  already  described.  There  are,  however,  considerable 
practical  difficulties  attendant  upon  designing  a  springing  arrangement 
to  give  the  correct  static  tension,  to  center  properly,  and  still  to  be 
much  “softer"  dynamically  than  that  which  is  used  in  the  experimental 
system.  Second,  may  be  eliminated  by  making 

Z„  »  (y  +  ivy'1  ,  (71) 

leaving  a  pure  traveling  wave  on  the  wire.  Reference  to  Cq.  (68)  and  con¬ 
sideration  of  the  magnitudes  involved  show  that  this  may  be  achieved  if 

a  M  Z  G  (72a) 

%  %  SG  (y  +  G)  .  (72b) 
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The  result  is  independent  of  frequency,  but  the  values  of  spring  and 
damping  constants  required  in  the  suspension  in  order  that  Eq.  (72a, b)  be 
satisfied  depend  upon  airspeed  and  altitude.  Such  an  arrangement  is 
therefore  of  rather  doubtful  utility.  It  should  be  remembered  also  that 
the  method  of  suspension  described  above  fails  entirely  to  eliminate  the 
longitudinal  motion. 


FIG.  12 

LONGITUDINAL  AMPLITUDE  AT  LOWER  END  OF  CABLE  WITH  FIXED  INPUT  ACCELERATION 
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FIG.  13 

TINS  I  ON  FLUCTUATION  IN  CABLE  WITH  FIXED  INPUT  AMPLITUDE 


FIG.  14 

TENSION  FLUCTUATION  IN  CABLE  WITH  FIXEO  INPUT  ACCELERATION 


V  SUMURY  AID  CONCLUSIONS 


In  summary  it  can  be  said  that  the  aerodynamic  characteristics  of 
trailing  vires  at  both  subsonic  and  supersonic  speeds  are  well  understood. 

In  particular,  the  factors  leading  to  cable  failure  are  knovn  and  have  been 
subjected  to  mathematical  analysis.  Investigation  of  various  means  for 
control  of  the  undesirable  phenomena  indicates  a  strong  possibility  that 
trai ling-wire  antennas  may  be  practical  even  for  supersonic  aircraft. 

The  static  shape  of  a  trailing  cable  can  be  described  very  closely 
for  all  speed- altitude  combinations  by  a  single  dimensionless  expression. 

The  mechanical  waves  which  can  exist  on  the  wire  have  also  been  described 
in  dimensionless  form.  The  theory  developed  in  this  report  expresses  the 
motion  in  terms  of  an  amplified  wave  traveling  downwind,  a  damped  wave 
traveling  upwind,  and  a  standing  wave  resulting  from  wire  curvature  and 
existing  primarily  in  regions  where  the  curvature  is  large.  On  this  same 
highly  curved  portion  of  the  cable  there  is  a  tendency  for  the  transverse 
oscillations  to  give  rise  to  longitudinal  oscillations.  If  the  longitudinal 
motion  is  restrained  in  any  way,  first-order  fluctuations  in  tension 
result. 

The  lower  end  of  a  trailing  antenna  can  be  stabilised  through  the  use 
of  a  heavy,  aerodynamic&lly  stable  end-weight  of  low  drag.  Under  these 
conditions,  the  most  important  vibratory  inputs  are  those  at  or  near  the 
upper  end  of  the  wire,  and  the  crucial  factor  determining  system  instability 
is  the  aerodynamic  amplification  experienced  by  the  downwind  wave.  The 
high  value  of  cable  tension  limits  the  amplification  to  a  relatively  low 
figure,  but  cannot  eliminate  it  entirely  at  supersonic  speeds  without 
exceeding  the  breaking  stress  of  the  cable. 

The  theory  described  in  this  report  was  applied  to  the  design  of  an 
experimental  system  which  has  been  successfully  flown  at  speeds  as  high  as 
Mach  1.24.  The  theoretical  amplification  in  the  worst  flight-test  condition 
was  about  8,  and  the  tension  fluctuations  were  calculated  to  be  of  the  order 
of  20  percent  for  liberally  assumed  vibratory  inputs.  The  tests  showed 
clearly  that  this  degree  of  aerodynamic  amplification  is  tolerable,  as  the 
observed  tension  fluctuations  attributable  to  the  phenomena  discussed  in 
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this  report  were  relatively  slight.  Tension  variations  resulting  from 
other  causes  were  observed,  but  none  was  of  sufficient  magnitude  to  bring 
about  failure  of  the  antenna  in  flight.  A  practical  antenna  for  use  on 
Mach  2  to  Mach  3  aircraft  would  have  maximum  amplification  of  the  same 
order  of  magnitude  aa  the  experimental  system.  Although  the  indications 
are  favorable,  the  question  of  the  stability  of  such  a  higher-speed 
installation  can  be  finally  settled  only  through  appropriate  flight  testing. 
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EQUATIONS  OF  THE  CABLE  SBAPE  IN  TERMS  OF  THE 
DIIENSIONLESS  VARIABLES  £,  V>  AND  <f> 
FOR  ZERO-DRAG  END-WEIGHT 


VARIABLE 

IN  TERMS  OF  £ 

IN  TERMS  OF  ff 

IN  TERMS  Or  <r 

IN  TERMS  OF  0 

— 

in 

WEES  M 

esc  0  -  1 

V 

cosh”1  (1  A-  £) 

... 

•ink"1  c 

In  cot  (0/2) 

<T 

•inh  rj 

... 

cot  0 

0 

esc-1  (1  ♦  £) 

cot"*  (sink  t)) 

COt~*  cr 

... 

S3 
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APPEND  I X)  B 


VATE  NOTION  IN  IK  CABLE  AT  ZENO  ANGLE  OF  ATTACK 

1.  THE  NON-LINEAR  EQUATION  OF  MOTION 

The  theory  given  in  Sec.  Ill  for  wave  Motion  on  a  wire  i Merged  in 
an  airstreaa  is  not  applicable  at  very  snail  angles  between  wind  and  wire. 
The  purpose  of  this  Appendix  is  to  develop  an  approxinate  theory  which  is 
valid  when  the  wind  direction  is  exactly  parallel  to  the  cable.  Before 
proceeding,  however,  it  is  convenient  to  describe  a  Mathematical  develop¬ 
ment  ihich  will  be  needed. 

Consider  the  function 

y(x)  *  |ain  x|  sin  x  (B-l) 

shown  in  Fig.  B-l.  This  function  can  be  expanded  in  a  Fourier  series  as 
follows; 

OD 

y(x)  =  £  ak  sin  kx  (B-2) 

*‘l 
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Only  odd  bornonics  are  present,  and  their  coefficients  are  given  by 


8 

*  k  (**  -  4)  IT 

Nuaerically, 

8  8  8 
°  1  , 3 ,  5 .  .  .  “  3?r  ■  lg7r  •  105?7 


(B-3) 


(B-4) 


The  coefficients  of  the  higher  harmonics  decrease  approximately  as  the 
inverse  cube  of  k. 

Nov,  for  slope  angles  <p  in  the  neighborhood  of  <t> 0  *  0,  the  second 
term  of  Eq.  (22)  cannot  be  neglected  in  comparison  with  the  first,  and 
the  approximation  of  Eq.  (23)  does  not  hold.  The  angle  4> ,  however,  may 
be  approximated  as  follows: 

.  Bn 

<t>  ■  sin  <p  *  tan  <p  «  — —  ( B-  5 ) 

OS 


The  aerodynamic  force  on  an  element  of  cable  is  proportional  in  magnitude 
to  K*  and  has  the  same  algebraic  sign  as  V  .  It  may  be  written  in  the 
following  form: 


where 


dF  *  f{s,t)  T  ds 


/(*,t) 


(B-6) 


(B  7) 


Substitution  into  Eq.  (24)  gives  for  the  equation  of  transverse  motion 


B*n 

Bs1 


JL 

£/* 


B2n 

Bt* 


+  /<«.*> 


0 


(B-8) 


To  obtain  an  approximate  solution  of  (B-8)  in  a  particular  case  where 
damping  is  known  to  be  small,  an  undamped  solution  is  used  to  evaluate 
VM,  and  thence  the  small  term  f(stt).  The  damped  solution  is  then  computed 
from  the  resulting  equation. 
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2.  TRAVELING  WAVES 


Suppose  that  the  solution  desired  is  that  of  a  traveling  wave  in  the 
downwind  direction,  excited  sinusoidally  at  the  origin.  The  dependence 
on  tine  nust  be  harnonic,  and  a  solution  of  the  forn 


is  assumed,  where  u(s) 


n  •  u(s)  sin  o)  |t  +  —  j 

is  real  and  slowly- varying. 
^  n  Bn  1  Bn 

IT  m  Ts  ~70~b7 


Then 


(B-9) 


(B-10) 


or,  if  the  tern  in  the  derivative  of  u  is  neglected, 


Also , 


a hi 

*  \ 

V 

n 

“  (1  -  y)  sin 

*  ( 

t  +H 
(/» 

+  — 

2 

(B-  11) 


/(*,<) 


C 


(l  -  y)2  y 


(B- 12) 


which,  by  Eqs.  (B-2)  and  (B-3),  is  seen  to  consist  of  a  Fourier  series 
whose  fundamental  leads  the  displacenent  by  90  degrees. 

With  this  linearization,  Eq.  (B-8)  can  now  be  treated  more  con¬ 
veniently  with  complex  numbers.  Since  the  complete  solution  nust  also  be 
periodic  in  time,  although  higher  harmonics  may  be  present,  it  will  be 
taken  to  be  of  the  form 


n 


£ 
o  rf  4 


«*(») 


eika(  t+t/U) 


(B- 13) 


The  boundary  condition  at  the  origin  is 


"(0, t)  -  ng  e 


(B- 14) 
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•nd  therefore 


«»(0)  -  n# 

(*  - 

1) 

(B-  15a) 

«*(0)  -  0 

<*  f 

1) 

(B- 15b) 

Equation  (B-8)  can  now  be  written  as 
B*n  l  3*/»  imu  \i 


3*J  (j*  3t2 


y  a 

1  V  ; 


~  y) 


2  • 


Y  eiia ><•*•/») 


*»i 

odd 


3*"  1  3»n 
3s*  Us  dt* 


■  ‘  Rr) ’  u~rn'  i  .. 

'  '  *  *  J 

•  44 


(B- 16) 

e  iia(  t*  $/ v) 

(B- 17) 


**  are  given  by  Eq.  (B-3).  When  ( B- 13)  it  substituted  in  (B-17) 
the  resulting  equation  must  be  satisfied  component  by  component,  and 
have 


we 


d2 


+  i  2k  — 


a>  duk 


d*2  '  '  "  V  ^7  •  i  7  n  -r)2  fl- 


am , 


(B- 18) 


The  aaplitudes  ut  will  be  assuaed  to  be  aoffici 

d*u 


'*  o>  rf“* 

77  «  2* - 

«*»  *  1/  <fs 


iciently  slowly- varying  that 

(B- 19) 


Then  Eq.  (B-18)  can  be  rewritten  to  give  the  following  set  of  equations 


du 


— i  .  (-pXo+n  yd  ~  r)2 

ds  -  *  — 


2*0/ 


(B- 20) 


(0 


•  u 


1 


In  pmicvltr,  for  k  •  lr  Eq.  (B-20)  b«co»e* 


dm 


i  4**1  -  y)J  * 

rfr  3wCV  1 


(B-21) 


which  can  he  iiUgnted  directly  to  give 


»,(*)  r  4rt\  -  y) 

l  + - 

3,-nCV 


^ .  f, 

"•  L 


—  VV^  A  1 


(B-22) 


If  the  foregoing  i*  eehetiteted  in  Eq.  (B-20),  the  aaplitedee  for  k  /  1 
c  an  he  fonad: 


»,(*) 


.  (_!)«( t-I) 


*»(*» 


4^(  1  -  y)  *  n 


yrrCU 


(B-23) 


Ia  diaeaaioaleaa  oaita  the  coaplete  aolatioa  can  be  written: 

■*]"* }  f  ,-^rrj, 

'  /  .Ji 


£i*Wr«r) 

(B-24) 


Equation  (B-24)  ibovs  that  the  fundanental  component  of  the  ware  ia 
amplified,  although  not  exponentially,  with  distance  along  the  wife,  while 
ana 11  amounts  of  higher  harnonica  are  created.  A  sinilar  investigation 
of  an  npwind  ware  given,  under  the  aane  conditions, 


Wcr,r) 


€  <v(r-«y) 


| 1  *  [‘  *  »  a  *  r)’  nvj  '  |  £  (-dXii-ii 


**(*»  -  4) 


£<»v(r-tr) 


tB-25) 

fron  which  it  can  be  seen  that  this  wave  ia  attenuated  with  distance. 

01 


i 
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3,  STANDING  WAVE 


Because  of  the  nonlinearity  of  Eq.  (B-8),  it  ia  not  possible  to 
superpose  the  traveling* ware  solutions  given  in  Eqs.  (B-24)  and  (B-25)  to 
obtain  a  standing-wave  solution.  Let  the  case  where  conplete  reflection 
at  the  origin  occurs  be  considered.  The  initial  solution  will  be  taken 
to  be 


n  *  nm  sin  —  sin  cot 


where  nB  is  real,  positive  and  constant. 

Using  Eq.  (B-10),  one  finds 

V*  nma>  f  ,  cos  9  .  „  cos  1 

K#  U  [CO*  U  +  7  *1D*  U  J 


(B-26) 


sin  (cot  -  \ p)  (B-27) 


where 


tan 


■if 

rtan  ~u 


(B-28) 


In  this  case  the  fundamental  conponent  of  f(s,t)  lags  the  displacement 
by  the  angle  i pt  and  we  can  write 


n‘‘,)  ■  T  pff  [co*’  ^  *y'  if] 


y(wt  t  «/») 

' (B-29) 


The  complete  solution  will  be  taken  to  be,  in  conplex  notation, 

”  "  *^i[6* 8in  (*t)  +  “*(s)  005  ' 


B-  30) 


If  the  system  is  excited  sinusoidally  at  the  point  »  »  S — that  is,  if 


n(S,t)  .  nfl  e4"* 


(B-31) 


the  boundary  conditions  become 


«4(0)  ■  0 


(B- 32a) 
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~  u,(S)  cot 


cjS 

kT 


(k  -  1) 

(k  f  1) 


qjS  coS 

”«  sin  “77  +  uils)  coa  ~7 7 


U  '  U 

With  the  use  of  Eq.  (B-29),  Eq.  (B-8)  can  be  written 

24.J.I4  .  .ip£l*  L..  22  , 


(B- 32b) 
(B-32c) 
(B-32d) 


,  •  j  w* 

3a*  0»  3t*  -  C\U  J  |C°8*T+r  810  7 


£  «*k  eiiat  e-***  . 

(B- 33) 


odd 


As  in  the  case  of  the  traveling  wave,  substitution  of  Eq.  (B~30)  into 
(B-33)  leads  to  a  set  of  equations  for  the  individual  components: 


d1uk 


1  °*  \ 

CO 

duk 

/ 

coa  (A  — 

-  2k 

-  sin 

A  — 

l  u  1 

V 

ds 

1  ul 

a*  /n«aM* r  a*  9  cm ]  ... 

•  (t)  [“*  t  *y,i"  tJ  e' 

Let  the  function  nfc(s)  be  defined  as  follows: 

du 


(B- 34) 


Vh($)  «  VkAS)  +  * 

Then  Eq.  (B~34)  can  be  written 


*  -  cos 

—  coa*  k—  .  (B-35) 

as  Xu* 


dv 

ds 


.  r  ,  a,.  ,  ,  ™ir 

r  ■  't(t|  [c“  n * r  *u  tJ [••• ** ■ 


i  sin  hp 


cos 


(COX  I 

*t) 

(B-36) 


Since  vk(rrU/2co)  must  be  zero,  the  real  and  imaginary  parts  of  vk  are 

t  f 


vir(«) 
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~C 


1  /n,<wi  f  r  ,  «■>*  *  .  * 1  1, 4>«»  .. 
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JirVllu 
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co*  T  +  r 


•  in3 


0)5 
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cos  U  ~  sin  fct/)  d* 


(B- 37b) 


f  _ 

Now,  from  Eq.  (B-32a)  ufc(0)  must  be  zero,  and  therefore  Eq.  (B-35)  can  be 
integrated  to  give 


uk(S)  « 


rs 

«  t  tos  \ 

vk(s)  sec*  (*~)  d* 

0 


(B-38) 


The  evaluation  of  u4(5)  in  the  general  case  requires  two  successive 
numerical  integrations  for  each  value  of  y  which  is  of  interest.  However, 
there  is  one  important  special  case  for  which  the  problem  is  greatly 
simplified.  Since  the  integrands  of  Eqs.  (B-37a)  and  (B-37b)  have  symmetry 
and  antisymmetry,  respectively,  about  the  point  s  ■  n(J/2o)t  the  corresponding 
integrals  have  the  reverse  properties.  Consequently,  when  the  cable  length 
5  is  equal  to  an  integral  number  of  half-wavelengths  the  real  part  of 
u*( 5)  vanishes  in  the  integration  of  Eq.  (B-38)*  The  imaginary  part  of 
t)k(s),  unlike  the  real  part,  can  be  found  analytically,  and  uk  for  any 
value  of  y  is  computed  by  means  of  a  single  numerical  integration. 

The  amplitude  of  harmonics  other  than  the  first  cannot  be  found  in 
this  special  case  by  the  method  under  discussion,  because  of  the  appearance 
of  cot  7 r  in  Eq,  (B-32c).  (In  principle,  the  inclusion  of  higher  order 
terms  in  Eq.  (B-33)  would  have  yielded  this  information. )  Results  for  the 
fundamental,  however,  are  obtainable  and  will  be  given  here.  Integration 
of  Eq.  (B-37b)  for  k  ■  1  gives 


v 


1  i 


SunsB 

9tt CU(  1  -  y* ) 


(B-39) 


•nd,  taking  into  account  the  fact  that  the  poaitive  3/2-power  muat  alwaya 
be  uaed,  Eq.  (B-38)  reduces  to 


C4 


c 


q(yY  ■ 


(B-40) 


The  index  p  gives  the  number  of  half- were  loops,  end  q(y)  is,  defined  by 


9(r) 


16 


rw  /  2 


9rr(l  -  y *) 


[(cos2  6  +  y2  sin2  6)***  “  y2]  sec2  0  cW 


(B-41) 


Numerical  values  of  q(y)  are  shown  graphically  in  Fig.  B-2. 

The  information  of  greatest  interest  from  the  point  of  view  of  design 
is  the  relationship  between  the  excitation  nQ  and  the  first-harmonic 
response  n m.  In  the  case  of  half-wave  resonance,  Eq.  (B-32d)  shows  that 
rt0  *  tti(^)»  ®nd  wc  can  write 

*  */P  9  O')  ’  JC  |n0|  .  (B-42) 

The  standing  ware  amplitude  is  seen  to  be  proportional  to  the  square  root 
of  the  input  amplitude. 


FIG.  B-2 

THE  FUNCTION  qly) 
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VARIATION  OF  WIRE  LENGTH  DUE  TO  STAPDING  WAVE 


If  the  first-order  wave  on  the  wire  is  known,  it  is  possible  to 
obtain  an  approximate  expression  for  the  variation  of  over-all  wire  length 
with  tine.  Let  ds  0  be  the  length  of  an  element  of  the  distorted  wire, 
n  its  displacement  normal  to  its  equilibrium  direction,  and  ds  its  pro¬ 
jection  on  that  same  direction.  Then 


(C-l) 


Note  that  n  here  represents  the  actual  instantaneous  value  of  the  dis¬ 
placement,  not  its  complex  value,  which  is  meaningful  only  in  the  first- 
order  theory.  For  moderate  amplitudes,  Eq.  (C-l)  can  be  expanded  as 
follows: 


(C-2) 


For  a  single  cycle  of  a  sinusoid  of  amplitude  a  and  wavelength  X  on  a 
straight  wire,  it  im  easy  to  show  by  integration  of  the  above  that 


(C-3) 


Now  the  instantaneous  transverse  displacement  along  the  curved  wire  is 
given  by  the  real  part  of  Eq.  (17).  If  the  time  origin  is  chosen  to  make 
Oj  real,  and  total  reflection  is  assumed,  the  result  is 


fle[n] 


(C-4) 
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This  equation  can  be  rewritten  as  follows: 


C^o 

fle[n]  *  a(<rn,t)  cos  |  +  &  I  - 


(C-5) 


where 


and 


•  «j  [4j2  +  4j2  -  2i4jl  cos  2cot]%  (C-6) 


V'(°'o>t)  *  tan-1 


*2  +  ^ 


4,  -  Ax 


tan  cot] 


(C-7) 


At  any  given  instant  of  time,  Eq.  (C-5)  represents  a  sinusoid  with  ampli¬ 
tude  and  phase  varying  along  the  cable.  For  waves  of  sufficiently  short 
wavelength,  this  variation  will  be  negligible  over  one  wavelength,  and  the 
contraction  of  the  wire  can  be  considered  as  the  sum  of  a  number  of  small 
contractions  in  which  the  local  length  decreases  by  the  factor  given  in 
Eq*  (C-3).  Then  it  is  possible  to  write 


(C-8) 


or 


Z  » 


[ A j2  +  A“2  -  2 A"1  cos  2<ot]  daQ 


.  (C-9) 


The  first  two  terms  in  the  integral,  which  can  be  evaluated  numerically 
in  a  particular  case  if  desired,  are  not  of  any  physical  interest,  as  they 
represent  only  a  static  change  in  the  equilibrium  length.  It  is  the  third, 
or  time- dependent ,  term  which  gives  the  amplitude  of  the  fluctuation. 

This  amplitude  is  easily  seen  to  be 
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VALIDITY  OF  EQUATION  (41)  AND  ANALYTIC  EXPRESSIONS 
FOR  THE  HOMOGENEOUS  SOLUTIONS 


The  purpose  of  this  Appendix  is  to  determine  the  range  of  conditions 
under  which  Eq,  (41)  is  a  good  approximation,  and  then  to  find  an  analytic 
representation  of  the  integral. 

Differentiation  of  Eq.  (41)  gives 


(D-l) 


The  condition  that  Eq.  (41)  be  a  good  solution  is  evidently  that 


f' 

2  vf** 


«  1 


(D-  2) 


or  say 
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10’1 


(D-3) 


The  quantity  on  the  left-hand  side  of  Eq.  (D-3)  always  has  a  maximum  for 
some  value  of  o\  If  this  maximum  is  less  than  10”1,  then  the  approximate 
solution  is  valid  for  all  values  of  cr.  Based  on  this  fact,  it  is  possible 
to  find  for  any  value  of  y  a  critical  value  of  v  below  which  the  approxi¬ 
mation  of  Eq.  (41)  fails.  The  resulting  curve  is  shown  in  Fig.  (D-l),  and 
will  be  referred  to  as  Criterion  A. 

The  evaluation  of  the  integral  appearing  in  Eq.  (41)  is  greatly 
simplified  if  [1  “  2i  (y/v)(  1  +  cr2)“^3^  can  be  expanded  in  powers  of 
2(y/v)(l  +  cr2)“^.  The  first  few  terms  in  the  expansion  are 

//(O’)  ■  1  -  i(y/v){  1  +  O’2)  ^  +  -^-(7/v)2  (1  +  a2)  1 

+  i(y/v)3  (1  +  O’2)  w  (y/v) 4  (1  +  O’2) 

£  8 
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FIG.  D-l 

CRITERIA  RELATING  TO  THE  APPROXIMATE  SOLUTION  OF  EQUATION  (39) 


Careful  testing  of  the  real  and  imaginary  parts  of  this  series  yields  the 
result  that  for  convergence  at  the  worst  point  (or  «  0)  one  must  have 

v  >  2y  (Criterion  B)  (D-5) 

This  criterion  may  be  seen  from  Fig.  (D-l)  to  be  less  restrictive  than 
Criterion  A  for  values  of  y  less  than  0.67,  and  more  restrictive  for 
7  >  0.67. 

The  imaginary  part  of  (D-4),  when  integrated,  leads  to  multiplicative 
factors  in  the  amplitude  of  the  wave,  while  the  real  part  leads  to  additive 
terms  in  the  phase.  To  ascertain  how  many  terms  of  Eq.  D-4  are  needed, 
the  integral  must  be  evaluated  term  by  term.  The  result  is 


(D-6) 


h(v) 


y%  _! 

-  tan  1  <y 

2v 


_5  ll 

16  vi 


a 

i  +  cr2 


+  tan*1 


.  (D-8b) 


Now  for  certain  ranges  of  y  and  v  it  will  be  possible  to  approximate 
the  two  functions  of  Eqs.  (D-8a,b)  by  unity  and  zero  respectively.  The 
deviatipns  will  be  worst  for  large  a,  and  in  fact  g(o*)  and  h(cr)  approach 
limiting  values  as  cr  becomes  infinite: 
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exp 
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(D-9b) 


The  factor  g(<r)  will  be  considered  approximately  unity  if  the  exponent  in 
Eq.  (D-9a)  is  less  than  0.1  or,  equivalently, 


v  >  2.24  y3/i  (Criterion  C)  .  ( D- 10) 

On  the  other  hand,  a  term  in  h(cr)  will  be  neglected  if  it  is  less  than 
0.1  radian.  This  condition,  applied  to  the  two  terms  of  (D-9b)  respec¬ 
tively,  yields 

v  >  15.7  y2  (Criterion  D)  ~  (D-ll) 

v  >  1.70  yA/i  (Criterion  E)  .  ( D- 12) 


Criteria  C  and  E  iare  seen  from  Fig.  D- 1  to  be  less  restrictive  than 
Criterion  A  for  values  of  y  which  are  of  physical  interest,  while 
Criterion  D  is  more  restrictive  only  for  y  >  0.24.  Thus,  we  have  for  all 
cases  of  interest  where  the  approximate  solution  is  at  all  valid 


g(a)  -  1 

(D- 13a) 

h(<r)  *  (y2/ 2v)  tan-1a  , 

(D-13b) 

and 

UjCa)  «  |<r  +  /l  +  <r2  exp  <  iv  la  +  ( y/v  )2  tan-1 

(D-I4a) 

U2(<r)  a  |<r  +  /l  +  <r2  exp^“iv|<T  +  -—  ( y/v) 2  tan~loj^ 

(D- 14b) 
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APPROXIMATE  METHOD  FOR  THE  COMPUTATION  OF  «,(<?) 
FOR  THE  EXPERIMENTAL  SYSTEM 


The  computation  of  u3(<r)  depends  upon  the  finding  of  a  particular 
integral  of  Eqs.  (56a, b).  If  the  inequality  of  Eq.  (57)  holds,  then  this 
system  of  equations  can  be  written  as  follows: 


*"  +  v2x  =  ( 1  +  a2  )~l* 

(E- la) 

y"  +  v*y  u  (  1  +  a2  )'1Ax 

(E-lb) 

Over  the  range  of  values  of  o  which  are  of  interest  in  the  current  problem — 
i.e.t  0  <  <r  ^  19,  the  right-hand  aide  of  (E-la)  can  be  represented  quite  well 
by  the  following  exponential  function: 

(1  +  o2)'l/l  rn  0.708  e'0  i22a  +  0.292 

+  0.440  o  g-1.8»«c. 

- 

(E-2) 

Figure  (E-l)  shows  the  original  function  and  the 
for  purposes  of  comparison.  It  is  seen  that  the 

approximation  of  Eq.  (E-2) 
approximation  is  very 

good  in  the  most  important  region,  near  the  lower  end  of  the  wire,  while 
it  deteriorates  somewhat  for  values  of  above  about  10. 

Uae  of  Eq.  (E-2)  in  Eq.  (E-la)  leads  to  the  following: 
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FIG.  E-l 

EXPONENTIAL  APPROXIMATION  FOR  + 

It  is  then  possible,  for  each  value  of  v,  to  obtain  graphically  an  approxi 
mate  expression  for  (1  +  <r2 )* l/i  x(a) „  For  the  case  at  hand,  an  expression 
of  the  form 


(1  +  o2yl*%(o)  -  (E-4) 

yields  sufficiently  good  results.  Substitution  of  Eq.  (E-4)  into 
Eq.  (E-lb)  leads  to  the  following  expression  for  y(<*): 


y(cr)  =  - — -  g-k  ( v  )<r  (E_5) 

[b{v)]  2  +  v1 

The  complete  solution  u3(<r)  is  then  easily  obtained  from  Eq.  (55)* 
Results  for  v  *  lt  2,  3  are  shown  in  Fig.  8  of  the  main  body  of  this 
report. 
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